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ABSTRACT 


Classical  (Langmuir)  probe  theory  is  reviewed  and  the  conditions 
under  which  it  is  applicable  to  the  re-entry  physics  problem  identified. 
Some  later  developments  in  probe  theory,  which  account  for  flow  velocity 
in  the  f ree - molecul ar  case,  are  presented. 

Measurements  of  current  densities  have  been  made  with  ion  probes 
mounted  flush  in  the  walls  of  a  pressure-driven  shock  tube;  the  results 
are  interpreted  using  a  one  -  dimens iona 1 ,  collision-dominated  probe  theory, 
allowing  for  the  spatial  nonuniformities  in  ion  density.  Correlation  of 
the  probe  theory  with  a  zero-order  gas  dynamic  boundary  layer  theory  is 
shown  to  lead  to  theoretical  current  estimates  in  fair  agreement  with 
the  experimentally  observed  values. 

Some  shock  ionization  precursor  results  in  air  are  also  presented. 
Finally,  the  results  of  some  calibrations  on  probes  that  were  flown  on 
the  RVIP  and  LORY  programs  are  presented. 
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I  INTRODUCTION 


The  determination  of  'the  electron  density  surrounding  a  re-entry 
vehicle  from  gas-dynamic  calculations  is  hampered  by  uncertainties  in  the 
rate  constants  of  the  ionization  and  de- ioni zation  processes  as  well  as 
by  the  effects  of  ablation  products.  However,  calculations  have  been 
made  for  ionization  in  air  (t.e.,  with  no  ablation  products);  and  it 
would  be  desirable  to  be  able  to  compare  these  calculations  against  flight- 
test  data. 

One  common  method  of  determining  electron  density  is  to  measure  the 
absorbtion  and  phase  shift  of  microwaves  that  propagate  through  the  plasma. 
However,  these  measurements  give  only  limited  information;  only  integrated 
values  for  absorption  and  phase  shift  can  be  obtained.  In  order  to  in¬ 
terpret  the  results  in  terms  of  electron  density  and  collision  frequency, 
it  is  necessary  to  estimate  the  spatial  distributions  of  these  parameters, 
although  the  spatial  distribution  is  itself  one  of  the  parameters  that 
one  would  like  to  obtain  from  measurements.  In  addition,  the  dynamic 
range  of  parameters  that  can  be  measured  by  such  methods  in  a  realistic 
flight  test  is  usually  small.  For  a  given  microwave  frequency,  it  is 
usually  possible  to  determine  the  electron  density  only  when  the  plasma 
frequency  is  of  the  order  of  the  microwave  frequency. 

An  alternative  technique  for  measuring  electron  density  during  flight 
tests  is  to  use  an  electrostatic  probe.  These  probes  can  be  designed  so 
that  measurements  can  be  made  with  a  spatial  resolution  of  less  than  1  mm. 
Thus,  the  use  of  a  probe  that  could  move  radially  outward  from  the  re¬ 
entry  vehicle  would  make  it  possible  to  determine  the  spatial  distribution 
of  electrons.  (One  could  also  use  several  probes  mounted  at  different 
radial  distances  from  the  surface.)  Electrostatic  probes  are  inherently 
capable  of  measuring  a  wide  range  of  electron  densities,  since  their  use 
involves  only  the  measurements  of  voltages  and  currents  that  change  rela¬ 
tively  slowly.  ' 

Classical  probe  theory,  as  developed  by  Langmuir,  is  predicated  upon 
a  number  of  conditions  that  will  not  hold  true  over  an  entire  trajectory.1* 


* 

References  ere  listed  nt  the  end  of  the  report. 
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It  is  necessary  to  determine  the  errors  involved  in  the  inferred  plasma 
parameters  and  to  develop  a  theory  that  will  enable  interpretation  of 
the  data  over  the  complete  trajectory. 

In  this  report,  the  classical  Langmuir  probe  theory  has  been  reviewed 
and  the  conditions  under  which  it  is  applicable  pointed  out.  Recent 
developments  in  probe  theory  that  account  for  flow  velocity  in  the  free- 
molecular  case  are  subsequently  presented  in  this  report. 

Preliminary  measurements  have  been  made  with  flush  ion  probes  mounted 
in  the  walls  of  a  pres  sure  -  dr i ven  shock  tuue.  In  this  position,  the  flush 
probes  are  exposed  to  the  large  ion-density,  temperature,  and  velocity 
gradients  existing  in  the  boundary  layer  in  the  vicinity  of  the  cold  wall. 
The  measured  current  densities  are  interpreted  using  a  one -dimens iona 1 , 
collision-dominated,  probe  theory,  allowing  for  the  spatial  nonuniformities 
in  ion  density.  Correlation  of  the  probe  theory  with  a  zero-order  gas 
dynamic  boundary  layer  theory  has  led  to  theoretical  current  estimates  in 
fair  agreement  with  the  experimentally  observed  values. 

Some  shock  ionization  precursor  results  in  air  are  also  presented. 

Finally,  the  results  of  some  calibrations  on  probes  that  were  flown  on 
the  Re-entry  Vehicle  Instrumentation  Program  (RVIP)  and  Low  Observable 
Re-entry  Vehicle  (LORV)  programs  are  given.  As  a  result  of  these  and 
many  other  measurements  using  probes,  it  is  concluded  that  measurements 
in  the  sa tura ted- ion- cur ren t  side  of  the  Langmuir  characteristic  are  con¬ 
siderably  more  accurate  than  measurements  on  the  electron-current  side 
near  plasma  potential. 
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II  ELECTROSTATIC  PROBE  THEORY 


A.  OPERATION  OF  PROBES  — CLASSICAL  THEORY 

An  electrode  or  probe  inserted  into  a  plasma  will  collect  charges 
until  it  has  sufficiently  high  potential  with  respect  to  the  plasma  to 
neutralize  the  net  charge  arriving  at  the  electrode,  i.e.,  the  electrode 
potential  is  such  as  to  equalize  the  electron  and  ion  flux.  If  the  elec¬ 
trode  potential  is  made  increasingly  negative  with  respect  to  a  second 
electrode,  fewer  and  fewer  electrons  will  reach  the  probe  until  only  pos¬ 
itive  ions  are  collected.  Conversely,  if  the  probe  potential  is  made  in¬ 
creasingly  positive  with  respect  to  the  second  electrode  (and  the  area  of 
the  second  electrode  is  greater  than  that  of  the  first  by  the  square  root 
of  the  mass  ratio  of  ions  to  electrons),  the  probe  will  collect  only  electrons. 

When  electrostatic  probes  are  used  to  collect  electrons,  the  satu¬ 
rated  electron  current  density,  J_,  is  related  to  the  electron  density, 
n_,  by  the  expression 


where 


n_e  v_ 

=  - 

4 


t/_  is  the  electron  random  velocity 

e  is  the  electron  charge 

k  is  Boltzmann’s  constant 
7L  is  the  electron  temperature 

is  the  electron  mass. 

Wh  en  ions  are  collected,  the  relationship  between  ion  current  density  and 
ion  density  (or  electron  density,  assuming  that  nj  *  n+)  has  been  shown 
by  Bohm  et  aZ.2  to  be  approximately 


J+  *  0.4 n+cv+ 


(1) 
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where 


is  the  velocity  with  which  ions  reach  the  sheath 
m+  is  the  ion  mass. 


From  the  expression  for  saturated  ion  current  density,  it  can  be  seen 
that,  if  the  ion  mass  and  electron  temperature  are  known  or  can  be  esti- 
mated,  a  measurement  of  the  ion  current  density  is  sufficient  to  deter¬ 
mine  the  electron  density.  Estimates  can  often  be  made  of  ion  species, 
and  hence  ion  mass;  in  addition,  the  electron  temperature  also  may  be 
estimated  from  gas-dynamic  considerations.  In  either  case,  errors  in 
these  quantities  will  not  be  very  serious,  since  they  both  enter  into 
the  expression  for  current  density  as  the  square  root.  For  example,  if 
the  temperature  at  a  certain  region  on  a  re-entry  vehicle  is  estimated 
at  3000°K  when  it  is  really  only  2500°K,  the  error  in  ion  density  is  only 
about  10  percent. 

Between  these  two  saturated  regions  is  a  region  that  may  be  used  to 
determine  the  electron  temperature. 

For  potentials  below  plasma  potential,  the  electrons  see  a  potential 
barrier  that  only  the  higher  energy  electrons  can  overcome  to  be  collected 
at  the  electrode.  If,  at  the  sheath  edge,  the  random  electron  current 
density  is  J_,  then  at  a  point  in  the  sheath  where  the  potential  is  V 
volts  below  the  plasma  potential,  the  electron  current  density  is 
J_exp(~eVWJ. 

The  net  current  density  is  the  algebraic  sum  of  the  ion  and  electron 
current  densities,  so  that  the  net  current  density  is 


Jnet  =  J_exp(-eK/feT_ )“J+ 


(2) 


Bringing  the  ion  current  density  to  the  left  side  of  the  equation  and 
taking  logarithms  of  both  sides,  we  obtain 

In ( J .  +  J  )  =  lnJ_  -  eV/kT _  .  (3) 

»  net 


Thus  T_  may  be  found  from  a  semi logari thmi c  plot  of  J4  +  J'  t  as  a 
function  of  V .  This  plot  will  be  a  straight  line  with  slope  e/feT.  for  a 
Maxwellian  energy  distribution. 
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Since  we  have  found  that  measurements  on  the  ion - satura tion- cur  rent 
side  are  more  accurate  than  measurements  on  the  el ec t ron - sa tu ra t i on  - 
current  side,  we  shall  restrict  the  following  discussion  to  the  ion  side. 

In  making  probe  measurements,  one  measures  current,  rather  than  cur¬ 
rent  density.  Therefore,  it  is  necessary  to  know  the  relation  between 
these  two  quantities.  From  a  consideration  of  the  formation  o  f  pi  asma 
sheaths,  the  sheath  thickness,  d  ,  can  be  approximated  by 

ds  -  \DV%  (4) 

where 


and  eQ  is  the  permittivity  of  free  space.  Equation  (4)  is  the  expression 
for  the  sheath  thickness  for  sp ace- ch arge- 1 imi ted  flow  in  a  planar  geo¬ 
metry..  This  equation  is  obtained  by  equating  the  current  density  given 
by  the  planar  space  charge  equation  to  the  random  current  density  of 
Eq.  (1).  The  approximation  of  a  planar  solution  for  sheath  thickness  is 
applicable  to  nonplanar  probes  only  when  the  sheath  thickness  is  small 
compared  to  the  probe  radius.  When  the  sheath  thickness  is  large,  it  can 
be  shown1  that  the  current  collected  by  a  probe  may  be  made  independent 
of  the  sheath  dimensions,  so  that  the  approximation  of  Eq.  (4)  is  no 
longer  required. 

The  overall  sheath  radius,  a,  is  given  by 

a  =  r  +  d  ,  (5) 

p  * 

where  rp  is  the  probe  radius.  If  the  probe  radius  is  larger  than  the 
sheath  thickness,  then  all  the  particles  that  pass  through  the  sheath 
edge  will  be  collected  by  the  probe  and  the  probe  area  will  be  approxi¬ 
mately  equal  to  the  sheath  area.  Under  these  conditions,  the  probe  cur¬ 
rent  for  a  cylindrical  probe  is 

I+  =  2nr pLJ +  (6) 
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and  for  a  spherical  probe 

J+  =  47rrJj+  ,  (7) 

where 

I+  is  the  probe  current 
L  is  the  probe  length. 

A  further  consequence  of  requiring  the  probe  sheath  to  be  of  the  same 
order  as  the  probe  radius  is  that  the  probe  current  will  be  relatively 
insensitive  to  the  probe  potential,  so  long  as  the  potential  is  suffi¬ 
ciently  negative  to  draw  saturation  current.  This  is  so  because  the  in¬ 
crease  in  probe  current  with  potential  is  due  to  an  increase  in  the  sheath 
thickness,  which  increases  the  area  through  which  current  passes.  How¬ 
ever,  since  the  sheath  area  is  approximately  equal  to  the  probe  area, 
there  will  be  no  significant  change  in  the  cur  rent- col 1 ec ting  area. 

In  order  to  ensure  this  condition,  it  is  necessary  that 

r  >>  d 
p  s 

and  that  the  probe  potential  be  such  that 


If  these  conditions  are  not  satisfied,  the  current  will  be  a  func¬ 
tion  of  voltage  and  will  no  longer  be  simply  equal  to  the  product  of  the 
random  current  density  and  the  physical  area.  Interpretation  of  probe 
current  in  terms  of  electron  densities  under  these  conditions  becomes 
more  complicated,  but  verified  theories  for  such  interpretations  are 
available.3-5  The  results  are  summarized  in  the  following. 

When  T)  >>  a/r  ,  where  a  is  the  sheath  radius  and  rp  is  the  probe 
radius,  the  current  is  independent  of  potential  but  is  a  function  of  a/r^: 

I+  =  ( a/r  ) 2nrpLJ+( for  a  cylindrical  probe) 


and 

1+  =  (a/ r  p)  2477T^J+  (for  a  spherical  probe). 

When  T)  <<  a/ r  and  rj  >  2 ,  the  current  is  a  function  of  potential  but 
independent  of  a/ r  . 

and 

J+  =  t?477T +  (for  a  spherical  probe). 

These  are  asymptotic  cases.  Hok*  has  calculated  the  dependence  of  ion 
current  on  ion  density  and  temperature  for  a  cylindrical  probe  for  the 
whole  range  of  r)  and  (a/rp).  For  the  spherical  probe,  the  intermediate 
range  of  rj  and  a/r^  can  be  approximated  using  the  solution  of  Allen,  et  aZ.5 

B.  CONDITIONS  FOR  APPLICABILITY  OF  CLASSICAL  PROBE  THEORY 

There  are  a  number  of  conditions  that  must  obtain  for  classical 
probe  theory  to  be  applicable. 

(1)  The  mean  free  path  must  be  much  larger  than  the  probe 
dimension.  This  is  necessary  to  ensure  that  the  probe 
does  not  disturb  the  plasma,  and  may  be  written  as 

rp  <  A.  (\  *  mean  free  path) 

(2)  The  sheath  thickness  must  be  much  less  than  a  mean  free 
path.  This  ensures  that  there  are  no  collisions  in  the 
sheath,  an  assumption  used  in  deriving  the  relations  be¬ 
tween  current  and  probe  potential.  This  maybe  written  as 

d  ~  a  -  r  <  k 

s  p 

(3)  There  is  no  net  flow  velocity;  this  may  be  written  as 


(4)  There  is  frozen  flow  across  the  sheath;  that  is,  there 
is  no  ionization  or  ionization  loss  — such  as  recombina¬ 
tion —  in  the  sheath.  This  may  be  expressed  as  a  current- 
continuity  relation 

V  J  =  0 


2irrpLJ+  (for  a  cylindrical  probe) 


7 


(5)  The  discussion  so  far  lias  concentrated  on  the  current 
collection  properties  at  only  one  electrode.  There 
must  be  another  electrode  in  contact  with  the  plasma 
to  complete  the  circuit  for  current  flow.  This  sec¬ 
ond  electrode  must  have  sufficient  area  for  current 
collection  that  the  sheath  around  it  does  not  limit 
current  flow  in  the  system.  This  consideration  of 
area  is  especially  important  in  making  the  tempera¬ 
ture  measurement  and  in  ensuring  that  conditions  are 
proper  for  measuring  the  random  electron  current  den¬ 
sity.  If  the  area  is  too  small,  a  saturation  current 
will  exist  for  both  ion  and  electron  currents;  however, 
the  electron  saturation  current  will  be  less  than  the 
full  random  electron  current  density.  The  area  ratio 
is  related  to  the  ratio  of  mobility  of  electrons  and 
ions;  in  air,  the  ratio  must  be  at  least  250. 

Area  Electrode  1 

-  >  250 

Area  Electrode  2 

If  the  el ec t rode- area  ratio  is  less  than  this  value, 
the  full  random  electron  current  density  will  not  be 
measured  at  any  potential.  The  electron  temperature 
can  be  inferred  from  the  cu rren t/ vo 1 tage  plot  but  the 
interpretation  will  be  somewhat  different  than  that 
outlined  above.  The  theory  of  probe  operations  at 
area  ratios  less  than  that  indicated  above  has  been 
worked  out  by  Johnson  and  Malter.6 

C.  ION  COLLECTION  WITH  FLOW  VELOCITY— FREE-MOLECULAR  CASE 

Several  theoretical  treatments7'9  have  investigated  the  case  of  ion 
collection  in  a  flowing  plasma  for  the  case  where  there  are  no  collisions 
in  the  sheath  and  where  the  probe  is  smaller  than  the  mean  free  path.  It 
is  also  usually  assumed  that  the  form  of  the  sheath  is  not  disturbed  by 
the  flow. 

The  theoretical  relation  between  probe  current  and  flow  velocity  is 
shown  in  Fig.  1  for  the  case  where  the  sheath  thickness  is  small  compared 
to  probe  radius.  Under  these  conditions,  a  useful  parameter  that  indi¬ 
cates  whether  the  flow  has  a  significant  effect  on  ion  current  collection 
is  the  ratio  of  flow  velocity,  v ^  to  thermal  velocity,  vf[vts  ( 2feT_/m+  Vz\  * 
This  ratio  is  approximately  0.8  times  the  Mach  number  for  air  before  dis- 
sociationl  When  the  Mach  number  is  much  less  than  unity,  the  effect  of 
flow  is  negligible.  When  the  Mach  number  is  greater  than  about  three, 
the  current  is  essentially  made  up  of  the  flow  of  ions  into  the  projected 
area  of  the  probe.  Thus,  for  a/ r  =1, 


FIG.  1  RATIO  OF  DIRECTED  TO  RANDOM  VELOCITY  AS  A  FUNCTION 
OF  THE  NORMALIZED  CURRENT  COLLECTED  BY  A  PROBE 
PERPENDICULAR  TO  THE  DIRECTED  VELOCITY 
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and 


I+  =  n+eVj2rpL  (for  cylindrical  probes  with  flow 
perpendicular  to  axis) 

I+  =  n  +  eVjTrr^  (for  spherical  probes). 

When  the  sheath  thickness  is  larger  than  the  probe  radius,  the  prob 
lem  is  more  complicated.  In  addition  to  the  Mach  number,  one  must  con* 
sider  the  relation  between  the  potential  field  around  the  probe,  which 
pulls  the  ions  toward  the  probe,  and  the  inertia  of  the  ions,  which  can 
carry  them  through  the  potential  field  without  being  collected.  This 
problem  was  treated  by  Smetana;9  his  results  are  shown  in  Fig.  2.  The 
actual  current  collected  is  obtained  by  multiplying  the  random  current 
collected  through  the  physical  area  of  the  probe  by  the  factor  F ,  given 
in  Fig.  2. 


FIG.  2  NONDIMENSIONAL  FLUX  DENSITY  OF  ATTRACTED  PARTICLES 
AS  A  FUNCTION  OF  DIMENSIONLESS  PROBE  POTENTIAL 
FOR  VARIOUS  SPEED  RATIOS,  S,  AT  a/r  «  100 
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When  S  (which  is  v  f/v  t)  is  equal  to  zero,  we  have  the  case  discussed 
previously  when  Tj  «  a/  r  .  The  current  increases  with  (17  +  1).  As  S 
becomes  greater  than  zero,  the  current  at  low  values  of  17  increases;  at 
high  values  of  7],  the  current  is  unchanged. 

This  behavior  can  be  explained  as  follows.  At  low  17  and  large  a/rpi 
the  potential  field  is  so  weak  that  many  of  the  ions  that  enter  the  sheath 
orbit  past  the  probe  and  leave  the  sheath  without  being  collected.  Under 
these  conditions,  the  current  collected  is  equal  to  the  random  current 
density  times  the  physical  area  of  the  probe.  Thus,  at  low  rj ,  when  the 
current  is  only  proportional  to  the  physical  area  of  the  probe  with  S  =  0, 
the  addition  of  flow  increases  the  current  collected,  just  as  it  would  in 
the  case  of  a/ r  ~  1.  At  large  values  of  S,  F  «  S  and  J+  »  nev^Apt  where 
Ap  is  the  projected  area. 

At  higher  values  of  Tj ,  with  S  =  0 ,  a  larger  percentage  of  the  ions 
that  enter  the  sheath  are  collected  by  the  probe.  From  simple  orbital 
calculations  one  finds  that,  with  a  given  potential  V  and  for  particles 
entering  with  initial  energy  VQ ,  the  radius  at  which  particles  will  just 
be  collected  is  r  =  r  (1  +  F/Frt)/2.  Therefore,  the  current  that  is  col- 

a  p  0 

lected  is  proportional  to  the  flux  entering  the  sheath  out  to  a  radius 
times  the  area  at  that  radius: 

I+  oc  nVj  x  rp  (1  +  V/V 0)‘A  . 

For  S  -  0,  F0  =  kT^/e  and  I+  is  given  by  the  equation  for  a/ r  >7}'  In 

V  p 

the  limit  of  7}'  »  1,  I+  ccV/2rp.  When  5  >  0,  the  flux  is  proportional  to 
SV ^  and  r  a  -  +  V/S2Fq)^.  Thus,  the  absorption  radius  is  decreased 

while  the  flux  entering  the  sheath  increases  so  that 

I+  oc  nSv'j  X  rp(  1  +  V/S2V0)‘A 

When  tj  »  S2  t  this  reduces  to  the  same  value  as  when  S  =  0. 

D.  ION  COLLECTION  BY  FLUSH  PROBES  IN  NONUNIFORM  PLASMAS 

In  order  to  avoid  perturbations  of  the  flowing  plasma  surrounding  a 
re-entry  vehicle,  it  may  be  necessary  to  place  probes  in  the  surface  of 
the  vehicle  (as  has  been  done  on  the  LORV  program),  rather  than  into  the 
flow  field.  Since  the  electron  density  varies  as  a  function  of  distance 
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from  the  vehicle  surface,  the  question  arises  as  to  what  current  the  probe 
will  measure. 

Ordinarily,  a  probe  samples  the  plasma  at  the  edge  of  the  sheath. 

For  the  classical  case,  where  the  mean  free  path  is  much  larger  than  the 
sheath  radius,  the  density  without  the  probe  present  cannot  change  signif¬ 
icantly  in  a  sheath  thickness.  However,  in  the  immediate  vicinity  of  the 
vehicle  surface,  the  electron  density  may  be  so  low  that  the  sheath  ex¬ 
tends  out  many  mean  free  paths  to  a  point  where  the  electron  density  has 
increased  significantly.  Under  these  conditions,  there  will  be  many  col¬ 
lisions  in  the  sheath. 

» 

The  following  is  a  first  attempt  at  a  on e- dimen sion al  theory  to 

handle  this  case.  s 

'i 

A  planar  geometry  is  assumed  with  the  probe  sufficiently  large  that 
the  physical  and  sheath  areas  are  approximately  equal;  most  of  the  cur¬ 
rent  drifts  into  the  sheath  from  a  plane  d  centimeters  from  the  surface. 
Assuming  that  there  is  neither  recombination  nor  ionization  in  the  sheath, 
an  equation  for  the  current  may  be  written  as: 

I  =  4n+  sevth/4  (8) 


where 


n  +  =  ion  density  at  the  sheath  edge 

v  fh  =  thermal  velocity  at  the  sheath  edge. 

If  these  assumptions  are  correct,  a  measurement  of  the  ion  current 
coupled  with  a  knowledge  (either  theoretical  or  from  another  measurement) 
of  the  thermal  velocity  at  the  sheath  edge  is  sufficient  to  determine  the 
ion  density  at  the  sheath  edge.  In  order  to  establish  the  location  of 
the  sheath  edge,  it  is  assumed  that  the  charge  density  redistributes  it¬ 
self  within  the  sheath  in  such  a  manner  that  the  current  flowing  across 
the  sheath  to  the  electrode  is  governed  by  the  same  relationships  as  for 
a  space- charge- limi  ted  diode  at  high  pressures.10  This  relationship  is 
given  by 

V2 

I  =  9/8  eoM+— 4  (9) 

d 3 
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who  re 


JjL +  =  ion  mobility 

K  =  potential  across  the  sheath 

d  -  distance  from  the  wall  to  the  sheath  edge 

A  =  area  of  the  electrode. 

Equation  (9)  may  be  solved  for  d$  in  terms  of  the  experimental  pa¬ 
rameters,!,  V  and  /x+ .  Equation  (9)  is  for  a  constant  mobility  through  the 
diode.  In  the  case  of  a  sheath  formed  near  the  vehicle  surface,  the  mo¬ 
bility  will  vary  through  the  sheath,  decreasing  toward  the  surface  by  a 
factor  generally  less  than  three.  The  diode  equations  have  been  solved 
with  a  mobility  decreasing  exponentially  from  the  cathode.  It  is  found 
that,  for  a  decrease  in  mobility  of  a  factor  of  three,  the  solution  has 
the  same  form  as  Eq,  (8)  with  /x+  equal  to  the  mobility  at  the  sheath  edge, 
but  with  the  constant  reduced  to  about  three  quarters  of  the  value  shown 
above.  Considering  the  order  of  the  other  approximations,  Eq.  (8)  will 
be  used  as  it  stands  with  equal  to  the  mobility  at  the  sheath  edge. 

If  these  assumptions  are  valid,  the  measurement  of  current  in  an 
electrode  flush  in  the  vehicle  surface  will  enable  one  to  determine  the 
ion  flux  at  a  distance  ds  from  the  surface.  It  is  interesting  to  note 
from  Eq.  (9)  that  the  sheath  edge  may  be  moved  by  varying  the  probe  po¬ 
tential,  thus  affording  the  possibility  of  measuring  the  electron  density 
profile  near  the  wall.  The  distance  ds  will  be  a  function  of  the  plasma 
parameters  as  well  as  the  potential.  For  reasonable  potentials  d!  will 
be  of  the  order  of  a  0.025  mm  (1  mil)  for  electron  densities  around  1013 
electron/ cc  and  increase  to  about  a  centimeter  for  densities  around  107 
electron/cc.  Thus,  the  degree  to  which  the  boundary  layer  near  the  ve¬ 
hicle  can  be  sampled  by  a  flush  probe  will  depend  strongly  upon  the  den¬ 
sity  levels  involved. 

Preliminary  results  using  this  theory  to  interpret  data  taken  in  a 
shock  tube  are  discussed  in  Sec.  III. 
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Ill  FLUSH  ION  PROBE  EXPERIMENTS 


To  date,  very  little  experimental  data  is  availabLe  regarding 
f lush- .noun ted  ion  probe  operation.  The  understanding  of  these  probes 
is  necessary  for  the  evaluation  of  the  data  from  the  RVIP  program.  The 
only  experiments  known  for  which  results  have  been  published  are  those 
of  Talbot11  and  Pollin,12  which  pertain  to  stagnation  point  flow  conditions 

The  experimental  determination  of  plasma  parameters  by  means  of  ion 
probes  mounted  flush  in  a  surface  bounding  the  plasma  necessarily 
restricts  the  spatial  extent  of  flush  probe  ion  collection  to  regions 
within  the  boundary  layer.  This  means  that,  for  cold  surfaces,  the 
flush  ion  probes  operate  in  spatial  regimes  where  large  ion  density  gra¬ 
dients,  thermal  gradients,  and  velocity  gradients  exist.  In  addition, 
because  of  the  relatively  low  ion  concentrations  in  the  immediate  vicinity 
of  the  wall,  the  probe  sheath  may  no  longer  be  free- molecular . 

A.  TEST  CONDITIONS 

To  obtain  preliminary  information  on  flush  ion  probe  operation  in 
high-velocity,  hi gh- tempe ra ture ,  free-stream  conditions,  an  experimental 
program  was  undertaken  at  SRI  in  the  12-inch  diameter,  pres s ure -dr i ven 
shock  tube.  In  this  exploratory  program,  ion  probes  of  various  geometries 
and  sizes  were  mounted  flush  in  the  shock  tube  wall.  In  this  position 
the  ion  probes  were  separated  from  the  free-stream  plasma  by  a  nonequi¬ 
librium,  nonsteady,  high-velocity  boundary  layer.  No  particular  care 
was  taken  to  remove  the  driver  impurities  from  the  shock  tube  walls  nor 
to  improve  the  tube  wall  roughness  (corresponding  to  a  standard,  12-inch 
diameter,  schedule  40,  aluminum  tube).  Untreated  room  air  was  used  as 
the  test  gas.  A  detailed  description  of  the  shock  tube  is  given  in 
Appendix  B  of  this  report. 

All  experiments  in  this  series  were  run  at  1  mmHg  initial  pressure, 
with  shock  velocities  ranging  from  2.8  to  6  mm/^sec,  The  corresponding 
range  of  theoretical  free-stream  ion  density  varied  from  1010  to 
2  x  1014ions/cc.  The  total  test  times  measured  from  the  arrival  of  the 
pressure  discontinuity  to  the  arrival  of  the  interface,  as  monitored  by 
photomultipliers  and  free-stream  ion  probes,  varied  from  zero  to  about  100/xsec. 
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B.  FLUSH  PROBE  MEASUREMENTS 

t 

Probe  geometries  and  the  associated  probe  circuitry  are  shown  in 
Figs.  3  and  4.  Probe  geometries  tested  included  a  9.5-mm  diameter 
circular  electrode  ^and  three  strip  electrodes  of  approximate  size 
6-3  x  1.25  mm  (HB-1),  7.5  x  1.25  mm  (HB-2)  and  6.3  mm  x  0.25  mm.  Since 
all  probes  were  connected  to  a  common  ground  through  the  oscilloscopes, 
they  were  operated  as  unequal  area  probes;  the  shock  tube  wall  served 
as  the  second,  c i rcu i t - c omp le t i ng  electrode,  Free-stream  electron  densi¬ 
ties  were  checked  in  a  separate  experimental  series  by  microwave  inter¬ 
ferometer  measurements.  Over  the  range  of  free-stream  conditions  covered, 
the  agreement  with  theoretical  equilibrium  values  corresponding  to  shock 
velocities  was  better  than  a  factor  of  two. 


FIG.  3  CONFIGURATION  OF  FLUSH  ION  PROBES 

(a)  9.5-mm  dia 

(b)  0.25  x  6.3  mm 
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0.23-  mm  DlA.,  12.5- mm  LONG 
PLATINUM  FREE  STREAM  PROBE 


0.025 -mm  DlA.,  WIRE  PROBE 


(d) 


TB- 3897-99 


FIG.  3  CONFIGURATION  OF  FLUSH  ION  PROBES 

(c)  HB-1 

(d)  HB-2 
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ION  PROBE 
ELECTRODE 


Typical  flush  ion  probe  responses  are  shown  in  Figs.  5  through  7. 
Figure  5  demonstrates  the  response  of  the  HB- 2  probe.  The  upper  trace 
corresponds  to  a  0.025-mm  (1-mil)  diameter  wire  probe  about  1.6  mm  above 
the  surface,  directly  over  the  flush  ion  probe.  The  wire  probe  axis  is 
perpendicular  to  the  flow.  The  lower  trace  indicates  the  response  of 
the  flush  ion  probe.  Note  that  the  flush  ion  probe  follows  the  wire 
probe  located  near  the  edge  of  the  boundary  layer  very  closely.  The 
current  level  of  the  flush  probe,  when  adjusted  for  probe  area  differences, 
is  about  a  factor  of  six  lower  than  the  current  level  of  the  wire  probe. 


STATION  579 
P,  *  I  mmHg 
Ut  *  4.4  mm/fiMC 
(n+)w  *  1*5  *  *0**  iont/cc 

0.5  v^m9IOOa 


0.02  y^mjooa 


U-MI7-K 


FIG.  5  TIME-RESOLVED  RESPONSE  OF  HB-2  ION  PROBE 
Upper  Trace:  Cylindrical  Probe  1.5  mm  Above  Wall  Surface 
Lower  Trace:  Flush  Probe 
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*  wmmmm 


Note  also  the  overshoot  in  the^flush  probe  current  at  the  arrival  of  the 
shock  front.  The  velocity  boundary  layer  thickness,  20  jJ* sec  after  shock 
arrival,  is  about  2  mm. 

In  Figs.  6  and  7,  characteristic  flush  probe  responses  to  widely 
differing  free-stream  conditions  are  shown.  Since  all  probes  indicate 
similar  response  characteristics,  only  the  HB- 1  probe  response  is  shown 
in  the  figures.  A  typical  probe  response  for  intermediate  and  high  free- 
stream  ionization  levels  is  shown  in  Fig.  6.  Note  the  characteristic 
current  overshoot  when  the  shock  arrives.  The  response  of  the  free-stream 

f 

ion  probe  of  0.25-mm  diameter  and  12.5-mm  length  is  indiated  in  the  same 
figure.  Figure  7  indicates  the  response  of  the  HB- 1  probe  to  free-stream 
ion  densities  about  one  and.  one-half  orders  of  magnitude  lower  ( llions/cc ) 
than  that  of  Fig.  6.  Both  the  flush  probe  and  the  free-stream  ion  probe 

— *j  — 20jA*ec 

STATION  579 
P\  =  I  mmHg 
Ut  =  4.25  mm/uMc 
(n+)«»  5  5  x  IO(f  iont/cc 

I  v/tm,IOU 


0.1  v/cm,  100a 

TA-SMT-tt 

FIG.  6  TYPICAL  TIME-RESOLVED  RESPONSE  OF  HB-1  ION  PROBE 

FOR  INTERMEDIATE  AND  HIGH  FREE-STREAM  ION  DENSITIES 
Upper  Trace:  Free-Stream  Probe 
Lower  Trace:  Flush  Probe 

show  the  slow  buildup  of  electron  density  behind  the  shock  and  the  overshoot 
of  the  flush  ion  probe  is  absent.  Note  that  in  both  figures,  corresponding 
to  high  and  low  free-stream  electron  densities,  the  f 1 ush- moun ted  ion  probes 
do  not  indicate  a  significant  increase  in  response  level  upon  arrival  of 
the  interface,  while  the  free-stream  probes  always  exhibit  a  substantially 
increased,  fluctuating  response.  Occasionally,  depending  on  free-stream 
turbulence,  the  flush  ion  probes  also  show  the  arrival  of  the  contact 
surface.  In  such  cases,  the  flush  ion  probes  follow  the  free-stream  probe 
response  closely. 
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I* — ZOfittc 


0.01  v/cm,300A 


STATION  579 
P,  =  I  mmHg 
Ut  *  3.36  mm/ustc 
(n+)w  s  1-06  x  I010  ionx/cc 


0.02  v/fcm,IOft 


FIG.  7  TYPICAL  TIME-RESOLVED  RESPONSE  OF  HB-1  ION  PROBE 
FOR  LOW  FREE-STREAM  ION  DENSITIES 

(a)  Upper  Trace:  Photomultiplier 
Lower  Trace:  Free-Stream  Probe 

(b)  Flush  Probe 
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Under  turbulent  free-stream  conditions,  am  generated  by  a  nonuniform 
shock  wave,  the  average  flush  proue  response  lefcel  is  substantially  in¬ 
creased  and  exhibits  rapid  fluctuations  (Fig.  Compared  to  a  laminar- 

type  prooe  response  for  identical  theoretical  free-stream  conditions,  the 
average  current  level  corresponding  to  turbulent  flow  is  about  a  factor 
of  two  higher. 


STATION  579 
P,  =  I  mmHg 
U«-  5.2 

(p+)m  s  8.5  x  10"  iont/cc 

2  v/cm,  lOA 


0.2  v/cm,  1000 


FIG.  8  RESPONSE  OF  HB-1  ION  PROBE  TO  TURBULENT  FREE-STREAM  CONDITIONS 
Upper  Trace:  Free-Stream  Probe 
Lower  Trace:  Flush  Probe 


The  current  overshoot  on  the  flush  probes  is  probably  related  to 
the  absence  of  a  developed  boundary  layer  immediately  behind  the  shock, 
the  translational  temperature  overshoot,  and  the  ion  density  incubation 
times.  The  overshoot  cannot  be  attributed  to  the  effects  of  probe 
capacity,  since  no  significant  correlation  between  probe  area  and  spike 
level  (or  between  probe  area  and  probe  time  response)  was  found.  Also, 
probes  covered  with  a  thin  coating  of  dielectric  (~10/x)  were  more  than 
one  order  of  magnitude  reduced  in  response  level,  and  the  transient 
response  was  significantly  shorter  than  observed  with  the  flush  ion  probes. 

C.  RESULTS 

The  experimental  probe  current  densities  have  been  plotted  agai  it 
the  free-stream  ion  density  (or  equivalently,  shock  velocity)  in  Fig.  9. 

The  theoretical  curve  appearing  in  this  figure  is  discussed  in  Part  D  of 
this  section.  Unless  the  estimated  relaxation  distances  behind  the  shock 
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FLUSH  PROBE  CURRENT  DENSITY - omp/m 


SHOCK  VELOCITY,  U, - mm /isec 


THEORETICAL  EQUILIBRIUM  FREE-STREAM  ION  DENSITY - ions/cc 


TS-JSST-60 


FIG.  9  EXPERIMENTAL  AND  THEORETICAL  FLUSH  ION  PROBE 

CURRENT  DENSITIES  vs.  FREE-STREAM  EQUILIBRIUM  ION  DENSITY 
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front  extended  further  than  10  cm  into  the  test  slug,*  the  data  were 
reduced  to  the  test  time  corresponding  to  10  cm  behind  the  pressure 
disc  ontinui ty . 

Interpreting  the  flush  ion  probe  data  on  the  basis  of  a  zero-order 
probe  theory,  described  in  Sec.  D,  it  was  found  that  the  flush  probes 
indicated  ion  density  levels  much  closer  to  free-stream  values  than  had 
initially  been  expected.  The  initial  estimates  were  based  on  the  assump¬ 
tion  that  the  cold  gas  near  the  wall  was  in  ionization  equilibrium.  At 
the  low- free- s tream  ion  densities  (1010  ions/cc),  the  ion  density  inferred 
from  flush  probe  data  was  a  factor  of  three  lower  than  the  corresponding 
free-stream  value.  At  the  higher  free-stream  ion  densities  (1014  ions/cc), 
the  flush  probe  data  indicated  values  about  an  order  of  magnitude  low. 

Presumably  because  of  poor  experimental  conditions  in  the  unsteady 
boundary  layer;  i.e.,  driver  impurities  deposited  on  the  shock  tube  wall 
and  turbulent  free-stream  conditions,  the  spread  in  experimental  data  is 
about  one  order  of  magnitude.  However,  when  observing  the  variation  of 
test  data  over  four  orders  of  magnitude  of  free-stream  ion  density,  from 
( 10 1 0  ions  /cc)  to  (1014  ions/cc),  it  becomes  evident  that,  for  the  present 
experimental  conditions,  the  flush  ion  probe  current  level  increases  with 
increasing  free-stream  ion  density  and  is  much  closer  to  free-stream  values 
than  considerations  of  ionization  equilibrium  at  the  wall  would  indicate. 

If  the  accumulated  data  from  all  probes  are  considered,  the  flush  probe 
current  varies  approximately  as 
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However,  a  better  fit  through  the  HB- 1  and  HB-2  probe  data  which  covered 
the  widest  range  of  free-stream  conditions,  would  yield 


J+  oc  n 


The  difference  in  the  power  dependence  is  caused  by  the  fact  that  the  larger 
(9.5-mm  diameter)  buttons  (used  largely  at  the  intermediate  free-stream  ion 
densities)  seemed  to  yield  higher  current  densities  than  Probes  HB- 1  and 


Here  the  test  slug  is  defined  as  the  region  of  shock  compressed  gss  between  the  shock  front  and  the 
interface* 
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HB-2  by  about  a  factor  of  two.  However,  no  definite  conclusion  can  be 
made  at  this  time,  since  this  difference  In  current  level  is  well  within 
the  experimental  data  scatter. 


D.  CORRELATION  WITH  THEORY 

To  check  the  experimental  results,  it  was  necessary  to  know  the 
variation  of  the  plasma  parameters  through  the  boundary  layer.  It  there¬ 
fore  became  necessary  to  investigate,  theoretically,  the  laminar  boundary 
layer  problem  for  the  flow  of  a  dissociated  gas  in  a  pipe.  To  obtain  a 
simplified  boundary  layer  solution,  a  coordinate  system  fixed  with  the 
shock  front  was  chosen  to  eliminate  the  time  dependence  of  the  boundary 
layer  flow. 

Since  laminar  flow  through  a  pipe  before  it  is  fully  developed,  can 
be  approximated  by  flow  over  a  flat  plate,  the  problem  of  a  uniform, 
compressible,  dissociated  gas  flow  over  flat  plate  with  a  sharp  leading 
edge  was  analyzed.  Analysis  of  a  flat  plate  boundary  layer  had  the  addi¬ 
tional  advantage  that  many  solutions  already  exist  providing  a  convenient 
check  on  the  computations. 

The  general  transformed  boundary  layer  equations  including  a  list 
of  symbols  is  given  in  Appendix  A.  These  equations  must  be  solved  by 
numerical  integration.  The  general  solution  for  the  shock  tube  case,  is 
currently  being  attempted. 

Initially  it  was  thought  better  to  investigate  the  behavior  of  a 
“  zero ’’-order  solution.  For  frozen  flow  in  the  boundary  layer,  a  cataly¬ 
tic  cold  wall,  i.e.,  z+w  =  0,  and  constant  coefficients,  the  ion  mass 
fraction  z+  is  given  by  Pr obs te in ' s 13  re  la tion : 

2+  =  /'(l?)  =  —  dO) 
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Since  the  ion  mass  fraction  is  defined  as 
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the  boundary  layer  ion  density  profile  is  given  by: 
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where  a  is  the  local  x  component  of  velocity,  and  /*  (r/)  is  relate^  to 
the  velocity  ratio  in  the  boundary  layer.  The  subscript  ®  denotes  ^'roc- 
stream  conditions.  Setting  the  Lewis  number,  Prandtl  number,  and  Schmidt 

number  equal  to  one  (these  parameters  are  defined  in  Appendix  A)  and 
assuming  constant  specific  heat  and  a  linearly  varying  viscosity  law,  the 
energy  equation  becomes  identical  to  the  energy  equation  for  undissociated 
air  as  solved  by  Chapman  and  Rubens te in 16 .  For  a  constant  wall  temperature 


and 
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The  functions  r{r])t  Y0(rj)}  r  (rj) ,  Y0(rj)  are  given  in  graphical  form  in 
Ref.  16  for  Pr  -  0.72.  These  functions  were  recalculated  for  Pr  =  1  for 
the  present  boundary- layer  solutions.  Typical  normalized  curves  for  the 
ion  concentration  profile,  velocity  profile,  temperature,  and  neutral 
density  profile  are  shown  in  Figs.  10  and  11.  Since  an  increase  in 
also  means  increasing  the  factor  /x^C,  while  y ^  and  are  approximately 
constant,  the  curves  of  Figs.  10  and  11  may  be  used  to  estimate  the 
boundary  layer  parameters  over  the  entire  range  of  test  conditions 
covered. 


Solving  Eq.  (8)  in  Sec.  II-D  for  ds  ,  one  obtains 
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where  J+  is  the  ion  current  density  and  the  subscript  s  denotes  variables 
evaluated  as  the  sheath  edge.  Since  the  numerical  value  for  the  mobility 
fj,s  depends  on  the  sheath  thickness,  Eq.  (15)  must  be  solved  by  iteration. 
However,  it  is  evident  from  the  1/3  power  in  the  exponent  that  ds  is  very 
insensitive  to  ion  mobility  and  ion  current  density.  Substituting  the 
experimental  flush  probe  current  densities  for  J+  and  estimating  the 
mobility  at  a  point  in  the  boundary  layer  corresponding  to  an  assumed 
sheath  dimension  from  the  surface,  an  improved  approximation  to  d$  can 
be  computed.  Once  the  solution  for  ds  converges  to,  say,  within  20  percent 
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NORMALIZED  TEMPERATURE  T/T, 


DISTANCE  OFF  PLATE  SURFACE,  y - mm 

TB-S857-97 


FIG.  10  NORMALIZED  TEMPERATURE  AND  MASS  DENSITY  PROFILES 

THROUGH  THE  LAMINAR  BOUNDARY  LAYER  ON  A  FLAT  PLATE 
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NORMALIZED  MASS  DENSITY  p/py 
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of  the  assumed  sheath  thickness,  the  approximate  sheath  thickness  has  been 
found.  The  sheath  dimension  for  the  range  of  experimental  conditions 
tested  is  shown  in  Fig.  12.  Note  that  the  sheath  dimension  becomes  very 
large  for  low  free- s tream  ion  densi ties  ex  tending  into  regions  where  the 
ion  density  at  the  sheath  edge  is  about  one-third  that  of  the  free-stream 
ion  density.  For  high  free-stream  ion  densities,  i.e.,  fast  shock  veloc¬ 
ities,  the  sheath  dimension  is  very  small,  extending  only  to  about  a  tenth 
of  the  free-stream  ion  density.  These  predictions  have  been  qualitatively 
verified  by  the  current  densities  obtained  in  the  flush  probe  experiments. 

For  a  given  sheath  dimension,  values  of  (n+)  and  ( V  ,)  can  be  obtained 

T  t  h 

s  s 

from  the  theoretical  boundary  layer  profiles  for  normalized  ion  density 
and  temperature.  Thus  theoretical  current  density  estimates  can  be  ob¬ 
tained  from  Eq.  ( 9 ) : 


<»♦>.<*, fc>.« 
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and  compared  to  experimentally  obtained  current  density  estimates.  The 
curve  for  theoretical  current  density  estimates  has  been  superimposed 
on  the  experimental  data  in  Fig.  9.  The  data  points  on  the  theoretical 
curve  correspond  to  sheath  dimensions  obtained  from  the  HB- 1  probe  data. 
Note  that  despite  all  assumptions  made  in  the  zero-order  boundary  layer 
theory  and  probe  theory,  the  agreement  between  theoretical  predictions 
and  experimental  data  is  surprisingly  good.  There  is,  however,  a  dis¬ 
agreement  in  the  power-law  dependence  between  the  experimental  and 
theoretical  current  relations.  The  difference  in  slope  is  probably  due 
to  the  mathematical  restrictions  on  transport  parameters,  frozen  flow 
condition,  etc.  A  priori,  there  is  no  physical  reason  why  the  free-stream 
ion  dens i ty - f lush  probe  current  density  relation  should  be  a  straight  line 
on  the  log-log  paper. 

As  expected,  agreement  is  best  at  the  low  shock  velocities,  where 
the  degree  of  ionization  is  Low,  while  at  the  high  velocities  the  approxi¬ 
mations  regarding  the  transport  parameters  and  the  frozen  flow  condition 
deteriorate. 
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EXPERIMENTAL  SHEATH  THICKNESS - mm 
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FIG.  12  EXPERIMENTALLY  DETERMINED  COLLISION-DOMINATED  SHEATH  THICKNESS 
FOR  FLUSH  ION  PROBE  vs.  FREE-STREAM  ION  DENSITY 

E.  SUMMARY  AND  CONCLUSIONS 

Over  the  range  of  free-stream  conditions  covered,  from  1010  to  about 
1014  ions/cc,  it  was  found  that  the  flush-mounted  ion  probes  sample  the 
ion  density  existing  at  the  probe  sheath  edge.  Interpreting  the  probe 
current  as  J+  =  n+evth/ 4,  at  the  low  free  stream  ion  density  levels,  the 
flush  mounted  ion  probes  measure  ion  density  about  a  factor  of  three  lower 
than  free  stream  values.  For  high  free  stream  ion  densities,  the  probes 
indicated  ion  densities  about  one  order  of  magnitude  reduced  from  free- 
stream  levels.  These  results  agree  well  with  theoretical  expectations 
based  on  boundary  layer  ion  density  profiles  and  sheath  dimensions  com¬ 
puted  by  means  of, the  collision  dominated,  space  charge  limited  diode 
equation.  In  addition,  considerations  of  probe  theory  indicate  another 
means  of  obtaining  boundary  layer  ion  density  profiles  from  flush  probe 
data.  Equation  (15)  shows  that  d*  has  a  strong  dependence  on  the  probe 
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potential.  Thus,  if  the  probe  sheath  can  be  varied  by  modulating  the 
probe  potential,  it  may  be  possiole  to  obtain  the  ion  density  profile 
in  the  boundary  layer. 

F.  FUTURE  WORK 

Since  the  flush  probe  environment  on  actual  re-entry  vehicles  is 
extremely  complex,  the  experimental  environment  has  been  simplified  in 
the  present  study  of  flush  probe  operation.  For  the  simplified  test 
conditions  the  variation  of  the  prone  parameters  through  the  boundary 

layer  becomes  calculable  and  the  effects  of  various  parameters  can  be 

\ 

isolated  and  studied  separately. . 

For  a  flat  plate  with  a  sharp  leading  edge,  suspended  in  a  uniform, 
high-velocity,  highly  dissociated  free  stream,  the  differential  boundary 
layer  equations  reduce  to  a  relatively  simple  form  and  numerical  integra¬ 
tion  is  not  too  difficult.  Accordingly,  the  next  phase  of  this  work  will 
be  primarily  concerned  with  flush  ion  probe  operation  on  flat  plates.  In 
addition,  whenever  possible,  the  theoretical  boundary  layer  calculations 
will  be  checked  experimentally,  by  measuring  the  boundary  layer  ion  density 
profile  thin-wire  probes. 

The  future  experimental  program,  as  it  is  presently  planned,  can  be 
described  as  follows: 

(1)  Determine  flush  probe  response  on  the  plate  over  a 
wider  range  of  free-stream  ion  densities, 

(2)  Determine  probe  response  associated  with  frozen  (as 
opposed  to  near - equi 1 ibrium )  flow.  These  parameters 
can  be  evaluated  in  the  shock  tube  by  changing  the 
initial  shock  tube  pressures. 

(3)  Determine  the  variation  in  probe  response  due  to 
changing  the  applied  battery  voltage.  This  also 
will  provide  a  check  on  the  probe  theory. 

(4)  Experimentally  determine  ion  density  profiles  in 
the  boundary  layer, 

(5)  Determine  boundary  layer  electron  temperatures. 

Currently,  experimental  work  on  Steps  (1),  (2),  and  (4)  is  in  pro¬ 

gress.  The  current  tests  are  being  conducted  at  0.1  mmHg  initial  shock 
tube  pressures,  covering  a  range  of  shock  Mach  numbers  from  Mg  =  16  to 
M$  =  24.  The  flat  plate  probe  configuration  currently  being  used  is 
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van  ous 


shown  in  Fig.  13*  Wires  about  0.025-mm  in  diameter  are  spaced  at 
distances  from  the  flat  plate  surface,  separated  laterally  to  minimize 
probe  interactions.  The  diameters  were  selected  so  that  the  probes  are 
in  a  free  molecular  regime  with  respect  to  the  neutral  gas.  To  avoid 
disturbing  influences  from  the  thicker  support  posts,  the  ends  of  the 
wires  are  insulated.  A  flush-mounted  probe  is  located  at  the  center  of 
the  plate.  All  probes  are  operated  as  unequal  area  probes. 


FIG.  13  EXPERIMENTAL  ARRANGEMENT  OF  CYLINDRICAL  ION  PROBES 
AND  FLUSH  PROBES  ON  A  FLAT  PLATE  TO  MEASURE 
ION  DENSITY  PROFILE  THROUGH  THE  BOUNDARY  LAYER 


Figure  14  indicates  the  nature  of  response  from  the  boundary  layer 
probes.  Note  that  small  perturbations  occurring  in  the  records  appear 
at  the  same  time  on  all  probes.  Also,  note  that  the  current  level 
fluctuations  are  small. 
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The  above  currents  have  been  reduced  to  ion  density  values  and  are 
shown  in  Fig.  15  superimposed  on  the  theoretical  boundary  layer  profiles. 
Note  that  the  agreement  between  theoretical  and  experimental  ion  density 
profile  is  good.  The  preliminary  result  should  be  viewed  with  caution, 
however,  since  it  represents  only  a  single  experimental  condition. 
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FIG. 
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15  REDUCED  EXPERIMENTAL  DATA  FOR  FIG.  14  SUPERIMPOSED 

ON  NORMALIZED  THEORETICAL  BOUNDARY  LAYER  ION  DENSITY  PROFILE 
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IV  PRECURSOR  MEASUREMENTS 


The  occurrence  of  precursor  ionization  phenomena  associated  with 
strong  shock  waves  in  air  and  in  the  noble  gases  has  been  reported  by 
several  investigators.1^16 

To  verify  these  reports,  experiments  have  been  performed  in  the  arc- 
driven  shock  tube  at  high  shock  velocities,  i.e.,  M$  =  15  to  22,  with 
room  air  and  at  0 . 1  mmHg  initial  shock  tube  pressures.  The  experiments 
used  a  free-stream  cylindrical  ion  probe  1.5  mm  in  diameter.  The  probe 
radius  was  of  the  order  of  a  free-stream  mean  free  path  ahead  of  the  shock. 

A  typical  ion  probe  response  is  shown  on  the  lower  trace  of  Fig.  16. 

It  is  interesting  to  note  that  the  rate  of  increase  in  signal  level  becomes 
very  large  a  short  distance  ahead  of  the  shock  front.  At  about  6  to 
10  /xsec  (7  cm)  ahead  of  the  pressure  discontinuity,  the  signal  is  driven 
off  scale. 

The  reduced  data,  adjusted  for  probe  sheath  dimensions,  are  presented 
in  Fig.  17.  The  bars  denote  the  extent  of  data  scatter.  For  reference, 
the  corresponding  equilibrium  slug  electron  densities  are  also  shown  on 
the  graph. 

To  verify  the  precursor  ionization  profiles  obtained  from  ion  probes 
and  to  discriminate  against  possible  probe  surface  photoelectric  effects, 

STATION  755 

P,  =  0.1  mmHg 

U,=  7.25  mm/fMC 

o.ooi  v/fcm,  10,000 a 

0.05  v/cm,  1000a 

FIG.  16  TIME-RESOLVED  MICROWAVE  INTERFEROMETER  AND  ION  PROBE  RESPONSES 
FROM  IONIZATION  PRECURSOR  MEASUREMENTS 
Upper  Trace:  3.3  Gc  Microwave  Interferometer 
Lower  Trace:  Free-Stream  Probe  (1.5-mm  diam  x  16  mm  long) 
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PRECURSOR  ELECTRON  DENSITY - el*ctrons/cc 
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TIME  BEFORE  SHOCK  FRONT  ARRIVED - u*«c 
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FIG.  17  MEASURED  PRECURSOR  IONIZATION  PROFILES 
IN  AIR  AHEAD  OF  THE  SHOCK  FRONT 
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a  few  experiments  were  repeated  using  a  33-Gc  microwave  interferometer. 

The  interferometer  actually  had  been  set  up  to  measure  test  slug  equi¬ 
librium  electron  densities  so  that  the  sensitivity  was  only  marginal  for 
the  precursor  measurements.  A  characteristic  precursor  microwave  inter¬ 
ferometer  response  is  shown  in  the  upper  trace  of  Fig.  16.  Note  that 
while  the  deflection  is  small,  it  is  observable  up  to  about  100  /x sec 
ahead  of  the  shock  wave.  To  reduce  the  probe  data,  the  temperature  of 
the  precursor  electrons,  has  been  assumed  equal  to  0.1  ev.  The  reduced 
microwave  data  have  been  superimposed  on  the  ion  probe  data  of  Fig.  17. 
Comparison  of  the  inferred  electron  densities  as  given  by  the  microwave 
system  with  those  inferred  from  ion  probe  measurements  show  agreement  to 
within  a  factor  of  two*  A  more  sensitive  interferometer  is  presently 
being  cons  true  ted  to  resol ve  the  electron  density  profile  in  the  precursor. 

The  effect  of  photoemission  from  the  probe  surface  on  the  probe  re¬ 
sponse  has  been  examined  by  enclosing  the  ion  probe  in  a  thin  dielectric 
sheath  (wall  thickness  =  0.3  mm).  Unless  the  pho t oemi s s i ve  properties 
of  the  dielectric  and  metallic  materials  are  identical,  the  current  re¬ 
sponse  of  the  dielectric-covered  ion  probe  to  the  precursor  should  change. 
If  the  current  is  not  produced  by  the  emission  of  photoelectrons  but  by  a 
plasma  of'charging  electron  density,  then  the  covered- probe  and  bare-probe 
current  responses  should  be  the  same  (the  capacity  of  the  dielectric  is 
very  large  compared  to  the  capacity  of  the  plasma  sheath;  hence,  there  is 
negligible  voltage  drop  across  it  as  the  electron  density  increases). 
Indeed,  the  bare  and  dielectric- covered  ion  probe  currents  were  found  to 
be  identical.  The  photoemissi ve  properties  of  the  dielectric  material 
(modified  polyvinylchloride)  are  unknown,  but  assuming  they  are  not  iden¬ 
tical  to  those  of  the  metal  conductor,  the  comparative  current  measurements 
indicate  that  photoelectron  current  is  negligible. 

The  present  results  indicate  electron  densities  of  the  order  of 
109  electrons/cc  about  10  cm  ahead  of  the  shock  front.  Verification  of 
the  precursor  ionization  is  important,  since  it  shows  that  radar  cross 
section  of  blunt  bodies  can  be  markedly  increased  at  VHF  frequencies. 

In  the  future,  it  is  planned  to  measure  the  electron  temperature  in 
the  precursor  by  a  Langmuir  probe  technique  and  to  measure  the  approximate 
distribution  of  nonequilibrium  shock  from  radiation  (assumed  to  be  causing 
the  precursor  ionization)  by  utilizing  the  photoelectric  effect  of  certain 
metal  surfaces17  and  a  field  of  view  perpendicular  to  the  direction  of 
motion. 
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V  RVIP  FLIGHT  TEST  CALIBRATIONS 


Recently,  Langmuir  probes  have  been  flown  on  the  Re-Entry  Vehicle 
Instrumentation  Program  (RVIP)  in  an  effort  to  measure  the  electron  con¬ 
centration  around  vehicles  re-entering  the  earth's  atmosphere.27  The  ini¬ 
tial  flight  results  indicate  that  the  inferred  electron  density  is  about 
an  order  of  magnitude  less  than  the  inferred  ion  density  and  that  the  in¬ 
ferred  electron  temperature  was  about  twice  the  value  expected.  In  an 
effort  to  resolve  some  of  the  flight  data  anomalies,  laboratory  calibra¬ 
tion  of  the  RVIP  probes  was  undertaken.  The  laboratory  measurements  have 
been  performed  in  an  electromagnetic  shock  tube  and  in  a  thermally  gen¬ 
erated  plasma.  The  environment  for  the  laboratory  measurements  was  such 
that  the  mean  free  path  was  much  greater  than  either  the  probe  radius  or 
the  sheath  thickness;  hence,  the  classical  theory  of  Langmuir  is  applicable, 

A.  MEASUREMENTS 

Three  aspects  of  probe  performance  are  to  be  checked  in  this  cali¬ 
bration  program.  It  is  desired  to  determine  the  applicability  of 
Langmuir  probe  theory  to  the  determination  of  electron  and  ion  density 
as  well  as  electron  temperature.  To  measure  the  electron  and  ion  satu¬ 
ration  current  characteristics,  the  electromagnetic  shock  tube  (described 
in  Appendix  C) ,  was  considered  most  suitable  because  of  its  ability  to 
reproduce  a  plasma  at  low  pressures  (long  mean  free  paths).  To  perform 
the  measurement  of  electron  temperature,  a  seeded  low-pressure  flame  was 
deemed  most  appropriate  because  of  the  long  test  times  available.  Details 
of  the  low-pressure  flame  facility  are  given  in  Appendix 

1.  General  Electric  Probe 

Figure  18  shows  one  of  four  sensing  tips  of  the  probe  designed  by 
General  Electric  mounted  on  a  support  column,  which  also  serves  as  the 
reference  electrode.  A  single  sensing  tip  and  the  column  serve  as  a 
Langmuir  probe.  The  area  of  the  reference  electrode  is  2,500  times  the 
area  of  the  sensing  tip.  The  sensing  tip  has  a  dimension  of  3 ( 10“ 2 )  cm 
and  an  area  of  4.2(10~7)  m2. 
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10>  (b) 

TA-MS7-IJ 

FIG.  18  RVIP  PROBE  GEOMETRIES  USED  IN  CALIBRATION 

(a)  General  Electric  Probe 

(b)  AVCO  Probe 

The  General  Electric  probe  is  designed  to  sample  the  hi gh- al t i tude 
(300,000  ft)  shock  formation  about  a  re-entering  vehicle.  As  a  conse¬ 
quence,  the  electron  density  regime  of  interest  is  between  107  and  109 
electrons/ cc.  In  this  electron  density  regime,  the  sheath  thickness  is 

much  greater  than  the  characteristic  probe  dimension.  The  details  of 
electrode  geometry  (i.e.,  whether  the  electrode  is  a  cube  or  a  sphere) 
are  thus  unimportant  and  the  electrode  is  assumed  to  be  spherical  for 
analytic  purposes.  (The  theory  for  the  collection  of  charge  by  a  spher¬ 
ical  probe  was  considered  in  Sec.  II-A. ) 

Figure  19  is  a  plot  of  ion  density  as  a  function  of  ion  saturation 
current  for  an  electron  temperature  of  5000°K  and  a  normalized  electrode 
voltage,  T)  =  eV/kT,  of  20.  The  solid  lines  of  Fig.  19  represent  the 
asymptotic  limits  of  current  for  probe  radius  much  greater  than  or  much 
less  than  ( T)  <  a/  r  )  the  sheath  radius.  The  open  circles  are  calculated 
points  using  the  theory  of  Allen  et  al.,5  and  the  dashed  lin e  i s  a  st rai gh t- 
line  approximation  joining  the  asymptotic  solutions  through  the  open 
circles. 

The  General  Electric  probe  was  placed  in  the  electromagnetic  shock 
tube  and  biased  for  ion  collection.  (Ion  density  can  be  inferred  from 
measurements  of  the  current  density  on  the  probe  as  the  shock  passes  the 
probe.)  The  ion  density  at  the  instant  of  shock  passage,  as  well  as 
400  fjise c  and  5  msec  after  the  passage  of  the  shock  front,  was  inferred 
from  curves  such  as  that  of  Fig.  19.  The  density  inferred  from  these 
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FIG.  19  ION  DENSITY  AS  A  FUNCTION  OF  ION  SATURATION  CURRENT 
FOR  THE  GENERAL  EL€CTRIC  PROBE 


measurements  was  checked  against  the  values  measured  by  a  cylindrical 
probe,  which  had  been  previously  calibrated  against  a  microwave  trans¬ 
mission  measurement.  Over  a  range  of  ion  density  from  10 7  to  109  ions/cc, 
the  probes  were  in  agreement  to  a  factor  of  two,  as  shown  in  Fig.  20. 

When  the  SRI  calibration  probe  and  the  RVIP  probe  were  biased  for 
the  measurement  electron  density  in  the  shock  tube,  the  inferred  electron 
density  was  an  order  of  magnitude  lower  than  the  inferred  ion  density. 

As  noted  previously,  the  error  in  electron  density  es timat es  i  s  f requen t ly 
observed  and  further  investigation  is  required. 

A  full  Langmuir  (current  voltage)  characteristic  of  the  General 
Electric  probe,  taken  in  a  seeded  flame,  is  shown  in  Fig.  21.  The  in¬ 
ferred  electron  temperature  from  the  slope  of  the  characteristic  is4000°K, 
a  number  that  is  a  factor  of  two  higher  than  the  estimated  temperature o f 
the  flame.  (A  similar  result  is  obtained  from  the  calibration  probe. ) 

For  this  case,  the  ratio  of  the  inferred  electron  and  ion  density  is 
within  a  factor  of  two  of  unity,  however  this  ratio  does  vary  with  elec¬ 
trode  voltage.  It  should  be  noted  that  the  plasma  environment  of  the 
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FIG.  21  LANGMUIR  CHARACTERISTIC  OF  GENERAL  ELECTRIC  PROBE 
IN  THE  FLAME  FACILITY 


seeded  flame  is  free  of  electric  and  magnetic  fields,  in  contrast  to  the 
environment  of  the  electromagnetic  shock  tube.  In  addition,  the  ioniza¬ 
tion  is  thermal  and  in  equilibrium  in  the  flame,  a  situation  that  does 
not  obtain  in  the  shock  tube. 

2.  AVCO  Flush -Mounted  Langmuir  Probe 

Figure  21(b)  shows  the  coaxial  geometry  of  the  Langmuir  probe  de¬ 
signed  by  AVCO  for  the  RVIP  program.  This  probe  is  designed  to  be 
mounted  flush  to  the  re-entry  vehicle  surface.  A  probe  having  the  same 
geometry  (but  slightly  different  dimensions)  of  that  of  Fig.  21(b)  was 
tested  in  the  electromagnetic  shock  tube. 

When  the  probe’s  inner  conductor  was  biased  negatively  for  ion  col¬ 
lection,  the  inferred  ion  density  agreed  within  a  factor  of  two  with  the 
measured  ion  density  inferred  from  a  calibration  probe,  provided  that  the 
collecting  area  of  the  flush  probe  is  taken  to  be  the  surface  area  of  the 
sheath.  The  sheath  thickness,  is  approximated  by 

ds  ~  X  £ 
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At  high  ion  densities,  when  the  sheath  thickness  is  smaller  than 
the  diameter  of  the  inner  conductor,  dt  the  area  lor  ion  collection  is 
just  the  physical  area  of  the  probe,  or  [ird^/A)*  At  low  ion  densities, 
when  d*  >  d ,  the  collecting  area  can  be  approximated  by  the  surface  area 
of  the  sheath,  or  (nd2s/2). 

When  the  inncer  conductor  was  biased  positively  for  electron  collec¬ 
tion,  the  inferred  electron  densities  were  too  low  compared  with  the  pre¬ 
vious  estimates  of  the  electron  density.  This  observation  is  in  agreement 
with  the  previous  measurements. 

One  difficulty  arising  from  the  geometry  of  the  AVCO  probe  is  its 
vulnerability  to  leakage  currents,  thereby  limiting  instrument  sensitivity. 
The  electrode  area  is  very  small.  At  low  altitude,  when  the  vehicle  sur¬ 
face  is  hot,  the  leakage  current  across  the  lossy  dielectric  can  quite 
possibly  exceed  the  current  collected  by  the  electrode,  particularly  at 
low  electron  density  levels. 

B.  CONCLUSIONS  AND  RECOMMENDATIONS 

On  the  basis  of  these  measurements,  it  is  concluded  that 

(1)  Determination  of  ion  density  from  the  RVIP  probes  is 
within  a  factor  of  two  of  the  actual  density  when  the 
classical  Langmuir  theory  is  applicable. 

(2)  Estimates  of  electron  density  from  the  RVIP  probes  are 
questionable.  This  conclusion  is  a  result  of  a  com¬ 
parison  of  inferred  electron  and  ion  densities  with 
the  electron  density  measured  by  microwave  techniques. 

In  these  measurements,  the  observed  electron  density 
can  be  as  much  as  an  order  of  magnitude  too  low. 

(3)  Estimates  of  the  electron  temperature  from  the  RVIP 
probe  data  are  no  better  than  a  factor  of  two  under 
ideal  conditions.  It  is  anticipated  the  accuracy  of 
the  temperature  measurement  will  suffer  when  the  elec¬ 
tron  density  estimates  differ  significantly  from  those 
of  the  ion  density.  In  any  case,  it  would  be  extremely 
desirable  to  record  the  entire  cu r ren t- vol t age  charac¬ 
teristic. 

An  examination  of  the  initial  probe  measurements  on  the  LORV  program 
has  led  to  the  following  recommendations: 

(1)  The  size  of  flush  probe  electrode  area  should  be  in¬ 
creased.  A  possible  configuration  requiring  a  minimum 
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of  system  changes  is  that  probe  being  used  in  a  joint 
program  involving  Sperry- AVCO- SRI .  (A  sketch  of  the 
configuration  of  this  probe  is  shown  in  Fig.  22.)  In¬ 
creasing  size  will 

(a)  simplify  interpretation  of  probe  currents; 

(b)  allow  measurements  at  high  altitudes; 

(c)  make  the  probe  less  vulnerable  to  leakage 
currents, 


FIG.  22  SPERRY  FLUSH-MOUNTED  PROBE 


(2)  Increase  the  voltage  swings  on  the  electrostatic  probe 
sweeping  circuitry  beyond  the  present  ±4  to  at  least 

±  10  volts.  The  probe  may  not  be  reaching  saturation 
electron  current. 

(3)  Operate  at  least  one  of  the  probes  with  a  constant 
voltage  bias.  A  suitable  bias  voltage  would  be  “10  v, 
so  that  probe  continuously  measures  ion  current.  This 
will  make  the  probe  less  susceptible  to  instrumenta¬ 
tion  difficulties. 

(4)  The  electrostatic  probe  electronics  should  be  checked 
out  in  a  plasma  environment,  to  check  for  any  pecu¬ 
liarities  introduced  into  the  results  when  the  instru¬ 
mentation  is  operated  with  a  nonlinear  load. 
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APPENDIX  A 


GAS-DYNAMIC  BOUNDARY  LAYER  THEORY 


1.  INTRODUCTION 

In  the  shock  tube,  the  flat  plate  boundary  layer  is  so  highly  cooled 
(constant  room  temperature  wall  conditions)  that*  (after  steady  state  is  • 
established)  the  presence  of  atoms  and  ions  in  the  boundary  layer  is  only 
due  to  diffusion  from  the  free  stream.  In  fact,  under  complete  tempera¬ 
ture  equilibrium  conditions,  the  charge  density  in  the  vicinity  of  the 
shock  tube  wall  would  always  be  equal  to  zero.  Considering  the  computa¬ 
tional  difficulties  associated  with  the  analysis  of  a  non-equilibrium 
multicomponent  gas  mixture  it  is  hoped  that,  for  the  shock  tube  conditions, 
the  simpler  boundary  layer  equations  using  a  binary  gas  and  average  gas 
mixture  properties  will  yield  useful  approximations  to  the  charge  density 
profiles  in  the  boundary  layer. 

The  following  general  assumptions  apply: 

(1)  The  gas  is  only  slightly  ionized. 

(2)  All  gas  species  can  be  considered  to  behave  as  perfect 
gases . 

(3)  The  flow  is  in  a  steady  state, 

(4)  Uniform,  parallel,  flow  over  a  flat  plate  is  considered, 
i.e.,  the  flow  is  two-dimensional,  and  all  free-stream 
gradients  in  the  x  direction  are  equal  to  zero. 

(5)  The  effect  of  a  boundary  layer  induced  oblique  shock 
over  the  plate  is  negligible. 

(6)  The  plate  temperature  remains  constant  at  the  room 
temperature  level,  i.e.,  the  plate  does  not  heat  up 
in  the  short  test  times  considered. 

(7)  Radiation  effects  are  ignored. 

Since  most  of  the  theory  used  here  has  been  developed  in  Refs.  18 
through  20,  only  major  steps  are  repeated  to  point  out  the  pertinent 
assumptions . 
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Cl 


The  following  symbols  are  used: 


=  speed  of  ™und 
=  Nondimens  ioLl 

=  STjt*  [(t,  + 


dens i ty- viscos i ty  product  =  ipp/p^p^) 

S)/{Tw  *  5)]  =  constant  of  proportionality  in 


the  absolute  viscosity  temperature  relation. 


c.  =  mass  fraction  of  component  i 
c+  =  mass  fraction  of  ions 

cp  .  =  specific  heat  per  unit  mass  at  constant  pressure  of 

component  i 


cp  *  average  specific  heat  per  unit  mass  at  constant  pressure 
for  the  gas  mixture,  defined  by  Eq,  (A-ll) 

Dx 2  =  binary  diffusion  coefficient  for  a  tom- molecu le  gas  mixture 

Z)q  -  ambipolar  diffusion  coefficient 

g  =  normalized  enthalpy  parameter  in  the  boundary  layer 


h  =  enthalpy  per  unit  mass  of  mixture 

h.  =  perfect  gas  enthalpy  per  unit  mass  of  component  i, 
defined  in  Eq.  (A-9) 

h®  =  heat  evolved  in  the  formation  of  component  i  at  0°K  per 
unit  mass 


hQ  =  total  enthalpy  defined  in  Eq.  (A-7) 
k  =  thermal  conductivity 

l  =  significant  dimension,  distance  from  plate  leading  edge 
L  =  Lewis  number  for  atom-molecule  mixture  =  (pD^  2.Cp  f /k ) 

M  =  Mach  number  =  U/a 


m  =  molecular  weight  of  mixture 

m+  38  molecular  weight  of  ion 

*-  molecular  weight  of  electrons 


n  =  neutral  particle  number  density 
n+  =  ion  number  density 

s  electron  number  density 
p  -  pressure 

Pr  s  Prandtl  number  =  \x/ kc 
Re  »  Reynolds  number  =  Ulp/fi 


45 


S  =  Schmidt  number  =  p/pDl2 

Sa  ~  Schmidt  number  based  on  ambipolar  diffusion  coefficient 
s  =  boundary  layer  parameter  defined  by  Eq.  (A- 15) 

T  =  absolute  temperature 
u  =  x  component  of  velocity 
v  =  y  component  of  velocity 

w t  =  net  mass  rate  of  formation  of  component  i  per  unit  volume 
and  time 

x  =  distance  along  plate 
y  *  distance  normal  to  surface 
2.  *  normalized  mass  fraction  of  specie  i 

zt  -  normalized  mass  fraction  ions 
T)  -  boundary  layer  parameter  defined  by  Eq.  (A-16) 
fx  =  absolute  viscosity 
y  =  ratio  of  specific  heat 
p  =  mass  density 

\p  =  stream  function  defined  by  Eq.  (A-12) 

SUBSCRIPTS 

oo  =  denotes  free  stream  conditions,  external  to  the  boundary  layer 
w  =  wall 

i  =  ith  component  of  the  mixture 
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BOUNDARY  LAYER  EQUATIONS : 


a.  Equation  of  Conservation  of  Specie 

If  it  is  assumed  for  a  slightly  ionized  gas,  that  charge  diffusion 
can  be  described  by  a  process  similar  to  chemical  specie  diffusion,  and 
if  in  addition,  the  assumption  of  a  binary  gas  is  made,  the  conservation 
equation  for  the  ion  mass  fraction  becomes: 


Bc+  3c  + 

pu  17  +  17 


•*  *  37  V10- 17 


3c. 


(A-l) 


where 


<< 


(A- 2) 


and,  since, 


n+n+  >:> 


(A- 3) 


n  ,m .  +  £  n  m 

+  +  ii 


nn  where  the  summation  is  over  all 

neutral  atoms  and  molecules 


Here  is  the  ambipolar  diffusion  coefficient,  and  the  two- componen t 
gas  is  considered  to  be  made  up  of  neutrals  and  charged  particles.  In¬ 
herent  in  the  above  equation  is  the  assumption  that  pressure  and 
temperature  diffusion  gradients  are  small  compared  to  specie  concen¬ 
tration  gradients. 

Summing  over  all  specie,  the  equation  for  conservation  of  mass 
becomes 


3(pu)  3 (pv) 

3*  3y 


(A- 4) 


b.  Equation  of  Motion 


The  momentum  equation  for  uniform  two-dimensional  flow  over  a  flat 
plate  becomes: 


pu 


(A- 5) 
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In  this  equation,  the  assumption  has  been  made  that  the  viscous  dissipa 
tion  gradients  in  the  x  direction,  within  the  boundary  layer,  and  the 
viscous  contributions  from  the  induced  y  component  of  velocity  are  of 
second  order. 


c.  Energy  Equation 

Since  it  was  assumed  that  the  gas  is  only  slightly  ionized,  the 
fractional  amount  of  energy  contained  in  the  ionization  mode  is  negligible 
and  therefore,  the  ionization  recombination  process  itself  cannot  influ¬ 
ence  the  temperature  or  enthalpy  profiles  to  a  large  extent.  This 
assumption  leads  to  the  uncoupling  of  the  specie  conservation  and  energy 
equation,  i.e,, the  variable  transport  coefficients  may  be  assumed  to  be 
independent  of  each  other  in  the  two  equations.  The  energy  equation  is 
therefore  written  on  the  basis  of  a  binary  gas  made  up  only  of  neutral 

atoms  and  molecules. 


In  terms  of  enthalpy: 


3hn 


*hn 


pu  —  +  fiv  — 

ox  oy 


—  i£ 
Bx  P 

r 


B(u2/2) 

3y 


where 

u2 

*0  -*  +  7 

h  =  2c.  h. 

i  '  1 

K  ■  JT  *  *! 


(A- 6) 

(A-  7) 
(A-8) 
(A-  9) 


In  terms  of  temperature: 

_  bt  _  bt 

pucp  -  +  pvc„  ~ 


Bu2 


+  2  w . ( h 0  “  h-  ) 


I 

i 


CPi^l2 


Be. 

I 

7~ 


B T 
^ y 


(A- 10 ) 
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and 


j* 

h 


Cf  frs?  ■  f  e‘Cf‘  • 


in  both  above  equations,  on  the  right-hand  side,  all  gradients  with  respect 
to  the  X  direction  have  been  dropped. 


d.  Transformation  of  the  Boundary  Layer  Equations 


The  above  mentioned  approximate  boundary  layer  equations  still  con¬ 
stitute  a  set  of  nonlinear  equations  that  cannot  be  solved  without  further 
approximations. 


In  order  to  transform  the  nonlinear  differential  equations  to  ordinary 
differential  equations,  which  are  separable,  the  technique  suggested  by 
Li  and  Nagamatsu15  is  used: 


1 p(*>y ) 


u 


V 


—  pudy  +  4j(x,  0  ) 

p  Jo 

Poo 

p 

p® 

p  B* 


(A- 12 ) 


IA-13) 


This  transformation  changes  the  equations  from  Cartesian  coordinates  to 
stream- function  coordinates  and  the  independent  variables  can  now  be 
separated  with  the  help  of  the  substitutions: 


<//(s,  T)) 

=  N(s)f(v ) 

u(S,  Tf) 

= 

c;  (s,  V) 

=  ci(s)zi  (17) 

h0(s,  Tj ) 

=  [/>0(s)g(7j)] 

e  tc. 

where  for  a  flat  plate, 


s  =  Moo /V*®* 


(A- 15 ) 


(A- 16) 
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With  the  assumption  of  boundary  layer  similarity,  the  s  dependence  of 
the  independent  variables  disappears,  so  that  the  boundary  layer  proper¬ 
ties  are  not  a  function  of  17  only. 


Here  the  similarity  principle  implies  that  the  boundary  layer 
parameter  profiles  are  similar  at  all  x  stations,  i.e.,  multiplication 
by  the  factor  vx/xRe^  reduces  the  profiles  from  all  x  stations  to  a 
single  profile  at  the  reference  station.  For  a  flat  plate  and  constant, 


dissociated,  free-stream  conditions  the  assumption  of 

a  similarity 

is 

valid.  The  new  transformation 

i 

relations  are: 

\  * 

=  *•/<  t,) 

* 

-  umf 

(t?) 

(A- 17) 

c  . 

i 

c^z 

i  0?) 

ho 

-  <*•> 

cg(v)  • 

The  variable  Tj  remains  the  same.  The 

boundary  layer  equations  for 

a 

flat  plate  with  uniform  external  flow 

finally  reduce  to: 

Equation  of  Motion: 

( cf“  )  ‘ 

+  ff" 

=  0 

(A- 18) 

Energy: 

OH  t/e  ■ 

Mr- 

l) 

i  (h  )<*>  1 

'  0 

(A- 19) 

Conservation  of  Specie: 

/  c  \  1 

2*i  + 

(r  ■*) 

II 

+ 

+ 

PPco^aoPaa 

(A- 20) 

whe  re 

c  =  ±  P  =  ^  . 

S  =  — 
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5  =  _  j 

f  P^12CP 

'ft**  k 

P°  12 
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The  boundary  conditions  for  a  cold,  catalytic  wall  are 


at  y  =  0 

*i(0)  =  2i. 

/( o)  =  /„ 

/'CO)  =  0 

g(0)  =  g„ 

g'(0)  =  g l 


Tj  =  0 


0 


CO  7J 

2  .  (00)  -  1 


/(<*>)  -  o 

/'  (°°)  -  1 
g(oo)  -  1 


% 


g'(°°)  -  0 


The  above  set  of  differential  equations,  in  the  general  case  of 
variable  coefficients,  still  does  not  yield  a  closed-form  solution,  but 
numerical  integration  is  possible  and  relatively  simple.  The  uncoupling 
of  the  energy  equation  and  the  equation  for  the  conservation  of  species 
reduces  the  computational  labor  required  considerably.  Aside  from  the 
general  assumptions  made  at  the  outset,  the  accuracy  of  the  above  bound¬ 
ary  layer  equations  for  the  present  application  is  probably  at  least  as 
good  as  the  accuracy  with  which  the  various  transport  coefficients  and 
the  rate  constants  are  known. 


E.  ZERO- ORDER  BOUNDARY  LAYER  APPROXIMATION 

If  one  sets  C  *  1,  P  ■  1,  L  =  1,  S  *  S  ■  1  in  the  transformed 
differential  boundary  layer  equations  [Eqs.  (A-18)  through  (A-20)]  and 
assumes  frozen  flow,  the  equations  reduce  to 

Women  turn: 

r  +  //"  =0  (A- 21) 


Specie  Concentration: 


Z+  +  fl' 


Energy  Equation: 


ft 

g 


+  /  g' 


0 


0 


(A- 22) 


(A- 23) 


51 


The  momentum  equation  is  now  essentially  in  the  form  of  t»he  well-known 
Blasius  boundary  layer  equation  for  incompressible  flow.  \For  the  Blasius 
differential  equation,  numerical  solutions  are  available  irv  tabulated21 
form. 

From  the  definitions  of  f  (77) ,  the  boundary  layer  velocity  profile 
is  given  as: 

—  =  f'iv)  (A- 24 ) 

“oo 

and  the  three  boundary  conditions  are: 

y  =  0  T)  =  0  y  =  0°,  77  =  00 

/( 0)  =  0 

/'( 0)  =  —  =  0  f(0 0)  .  —  *  1  . 

Uco  “co 

If  one  substitutes 


Z+  =  f'(v) 

g  =  f'iv) 


(A-25 ) 


into  the  specie  conservation  and  energy  equation,  one  reduces  these  equa¬ 
tions  to  the  Blasius  equation.  For  a  cold  catalytic  wall,  (Z+)  =  0  and 

=  /(tj)  ■  —  •  (A-  26) 

Uoo 

This  equation  is  Probstein's  second  integral.  It  indicates  that,  for  the 
above  assumptions,  the  ion  mass  fraction  in  the  boundary  layer  varies 
similarly  as  the  x  component  of  velocity.  By  definition, 


(A-27 ) 


From  this  relation,  the  ion  density  profile  in  the  boundary  layer  is 
obtained  as: 


n 


+ 


(A- 28) 
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A  similar  expression  can  be  found  for  the  enthalpy  profile;  however, 
for  this  application  one  is  more  interested  in  the  temperature  profile. 

If,  in  addition  to  the  previous  assumptions,  one  neglects  the  effect  of 
diffusion  of  atoms  on  the  temperature  profile  in  the  boundary  layer,  the 
differential  equation  for  the  normalized  temperature  profile  may  be 
expressed  as: 

P 

T*"  +  PrfT*‘  =  --j  (ya  -  1)A£  {f")2  ,  (A-29) 

where 

T *  =  T/T . 

This  equation  for  the  temperature  profile  has  been  solved  by  Chapman  and 
Rubesin18and  the  solution  is  given  by: 

Too  ~  1 

-  =  1  +  ■  m  “  'Mi  r(7?)  +  a0Y0(V )  ■  (A-30) 

rr  IT 

1  00 

The  functions  r(^)  and  YQ(ri)  are  given  in  Ref.  18  for  Pr  =  0.72.  These 
functions  have  been  recalculated  for  P r  =  1  for  the  present  purpose. 

The  appropriate  boundary  conditions  are: 

MO)  =  i  M°°)  =.° 

r(0)  =  /Pr  r(°°)  =  0 

The  temperature  profiles  have  been  computed  using  equilibrium  free  stream 
y  and  A/qj  values,  and  correspond  to  the  profiles  associated  with  frozen- 
flow  conditions,  The  error  in  estimated  probe  current  due  to  the  approx¬ 
imated  temperature  profile  is  probably  not  larger  than  a  factor  of  2 . 

The  effect  of  inaccurate  temperature  profiles  will  be  felt  more  strongly 
in  the  relationship  between  the  boundary  layer  parameter  7)  and  y,  and 
will  yield  thermal  boundary  layers  thinner  than  the  actual  boundary  layer. 

Near  the  cold  wall,  large  gradients  exist  in  both  the  temperature 
and  ion  concentration  profiles,  and  relatively  small  charges  of  y  cor¬ 
respond  to  large  changes  in  the  numerical  values  of  ion  concentration. 
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Following  Chapman  and  Rubesin,18  the  relation  between  the  boundary 
layer  transformation  parameter  7]  and  the  physical  distance  y  from  the 
flat  plate  surface  is  given  by: 


1  /  Pcduo>  f y  n 

T?  =  —  / -  —  dy 

2  J  Ml  CO  X^l  0  Pco 


( A- 3 1 ) 


now,  since 


( A- 32 ) 


r  t 

y  .  2/——  —  dn' 

\J  ucdPcq  od 
v  Jo 


( A-  3  3  ) 


assuming  that  m^/m  ~  1.  Substituting  for  T/T^  from  the  energy  equation 
and  integrating,  the  relation  between  y  and  7]  becomes 


/Moo  XCi  y  -  1 

2/  tj  +  — - —  W*  r  (17)  +  a0Y0(i7) 


(A-34) 


where  the  functions, 


rV 


r( 7)'  )dr) *  Y(t)) 


Yo  (V‘)dr,' 


(A-35 ) 


have  been  calculated  for  P =  1.  It  is  evident  that  an  underestimation 
of  local  temperature  will  yield  boundary  layers  that  are  too  thin*  This 
effect  is  felt  most  severely  near  the  wall,  where  the  integral  is  still 
small.  However,  underestimation  of  local  temperature  will  be  offset  to 
some  extent  by  the  assumption  of  a  constant  molecular  weight  across  the 
boundary  layer. 
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APPENDIX  B 

PRESSURE- DRIVEN  SHOCK  TUBE 


The  laboratory  simulation  of  plasma  phenomena  associated  with  plane¬ 
tary  re-entry  conditions,  such  as  voltage  breakdown  during  re-entry,  and 
communication  “blackout,”  operation  of  free-stream  and  f lush -mounted  ion 
probes  on  re-entry  vehicles,  etc.,  necessitate  the  existence  of  a  high- 
velocity  plasma.  In  addition,  studies  of  the  strong  shock  ionization 
precursor  effects  and  investigations  of  the  shock  front  nonequilibrium 
radiation  overshoot  require  the  production  of  very-high- veloci t^  shocks. 

For  short  test  times,  the  arc-driven  shock  tube  can  provide  strong  shocks 
and  equilibrium  plasmas  moving  with  velocities  up  to  about  40,000  ft/sec. 
One  general  advantage  of  shock-tube  testing  is  that,  if  equilibrium  test 
conditions  exist  in  the  test  slug,  the  gas  properties  are  a  function  only 
of  shock  velocity.  Thus,  an  accurate  measurement  of  shock  velocity  yields 
all  equilibrium  gas  properties. 

1.  DESCRIPTION  OF  THE  EQUIPMENT 

The  shock  tube  used  in  these  studies  is  of  the  arc-driven  type  and  is 
modeled  after  the  design  of  Camm  and  Rose.22  A  schematic  of  the  shock  tube 
is  presented  in  Fig.  B-l.  Figure  B- 2  shows  a  cross  section  through  the 
shock  tube  driver.  The  driver  is  1%-inches  in  diameter  and  12  inches  long. 
The  driver  is  connected  to  the  12-inch  diameter  driven  tube  by  means  of  a 
54  inch  long  circular  transition  section.  The  length  of  the  driven  tube 
has  recently  been  increased  from  25  to  36  feet. 

The  energy  source  is  a  100  kilojoule  capacitor  bank,  which  is  trig¬ 
gered  by  an  exploding  wire  technique. 

2.  INSTRUMENTATION 

Shock  velocity  is  measured  by  a  series  of  flush  probes  mounted  in  the 
tube  walls.  The  probe  current  develops  a  voltage  across  a  resistor,  This 
voltage  is  then  differentiated  to  give  a  spike  in  voltage  when  the  ioniza¬ 
tion  front  passes  the  probe.  The  voltage  spike  is  applied  to  the  vertical 
plates  of  a  raster  oscilloscope  (see  Fig.  B-3),  and  the  time  for  the 
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FIG.  B-l  SCHEMATIC  LAYOUT  OF  PRESSURE-DRIVEN,  ARC-HEATED  SHOCK  TUBE 


FIG.  B-2  CROSS  SECTION  OF  SHOCK  TUBE  DRIVER 


V* 


P, » I  mm 

(^i)av9  =  4.5fnn\/^l#C 
100  /itec/line 


TA-MS7-44 


FIG.  B-3  SHOCK  VELOCITY  MEASUREMENT  ON  RASTER  SCOPE 


shock  to  travel  a  known  distance  is  measured  from  the  oscilloscope  in 
terms  of  the  time  between  spikes.  The  accuracy  with  which  velocity  is 
measured  is  estimated  to  be  about  2  percent.  Typical  response  curves  of 
various  instruments  used  on  the  shock  tube  are  shown  in  Figs.  B-3 
through  B-7. 

Figure  B-4  shows  a  pressure  probe  trace  and  the  time  resolved  re¬ 
sponse  of  a  0,5  mm  diameter  platinum  free-stream  ion  probe.  Note  the 
indicated  flat  ion  density  level  in  the  test  gas  slug  and  the  increased 
probe  response  when  the  interface  arrives.  (The  increased  interface 


0.2  v/cm,  I  megohm 

STATION  456 
Pi  *  I  mmHfl 

U,  «  4.75  mm//AMC 
(n+)mm  *  3.7  g  10"  ione/cc 

0.2  v/cm,  1/20 

TA-MST-4S 

FIG,  B-4  COMPARISON  OF  FLUSH-MOUNTED  PRESSURE  PROBE  RESPONSE 
AND  0.25  x  6.3  mm  FREE-STREAM  ION  PROBE  RESPONSE 
Upper  Trace:  Pressure  Probe 
Lower  Trace:  Free-Stream  Ion  Probe 
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response  is  attributed  to  the  arrival  of  driver  impurities  having  lower 
ionization  potential  than  air.)  It  should  also  be  realized  that  the  ion 
probe  measures  only  the  local  ion  density  level,  so  that,  in  general,  one 
cannot  conclude  from  these  data  anything  in  regard  to  the  state  of  the 
gas  at  other  localities  in  the  test  slug.  On  the  other  hand,  this  state¬ 
ment  can  be  turned  around  to  say  that  for  experiments  involving  ion 
probes,  only  the  local  test  conditions  are  important.  This  fact  has  been 
utilized,  since  experimental  experience  has  indicated  that  test  times  in 
the  center  of  the  shock  tube  can  be  substantially  lower  than  test  times 
near  the  shock  tube  wall. 

Figures  B- 5  and  B-6  show  typical  photomultiplier  traces  at  1  mmHg 
and  0.1  mmHg  initial  shock  tube  pressure,  respectively.  In  Fig.  B-5  note 
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FIG.  B-5  TYPICAL  PHOTOMULTIPLIER  RESPONSE 

FOR  1  mmHg  INITIAL  SHOCK  TUBE  PRESSURE 
Upper  Trace:  Photomultiplier 

Lower  Trace:  Free-Stream  Ion  Probe,  0.25  mm  dia  x  6.3  mm  Long 


that  the  local  indicated  test  time  of  the  ion  probe,  located  outside  the 
photomultiplier  field  of  view,  is  about  twice  as  long  as  indicated  by  the 
photomultiplier.  The  test  time,  as  observed  at  a  given  shock  tube  station 
is  defined  ris  the  time  interval  from  the  arrival  of  the  shock  front  to  the 
arrival  of  the  interface.  It  is  also  interesting  to  observe  that  the  ion 
probe  response  reaches  a  steady  level  sooner  than  the  photomultiplier,  and 
then  remains  constant,  despite  the  fact  that  chemical  relaxation  is  still 
proceeding.  The  photomultiplier  field  of  view  is  directly  along  a  shock 
tube  diameter.  Collimation  is  achieved  by  two  0.15  mm  wide  slots  placed 
12  inches  apart.  The  aperture  slots  are  about  3.25  mm  long. 
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response  is  attributed  to  the  arrival  of  driver  impurities  having  lower- 
ionization  potential  than  air.)  It  should  also  be  realized  that  the  Lon 
probe  measures  only  the  local  ion  density  level,  so  that,  in  general,  one 
cannot  conclude  from  these  data  anything  in  regard  to  the  state  of  the 
gas  at  other  localities  in  the  test  slug.  On  the  other  hand,  this  state¬ 
ment  can  be  turned  around  to  say  that  for  experiments  involving  ion 
probes,  only  the  local  test  conditions  are  important.  This  fact  has  been 
utilized,  since  experimental  experience  has  indicated  that  test  times  in 
the  center  of  the  shock  tube  can  be  substantially  lower  than  test  times 
near  the  shock  tube  wall. 

\  Figures  B- 5  and  B- 6  show  typical  photomultiplier  traces  at  1  mmHg 

and  0.1  mmHg  initial  shock  tube  pressure,  respectively.  In  Fig.  B-5  note 
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FIG.  B-5  TYPICAL  PHOTOMULTIPLIER  RESPONSE 

FOR  1  mmHg  INITIAL  SHOCK  TUBE  PRESSURE 
Upper  Trace:  Photomultiplier 

Lower  Trace:  Free-Stream  Ion  Probe,  0.25  mm  dia  x  6.3  mm  Long 


that  the  local  indicated  test  time  of  the  ion  probe,  located  outside  the 
photomultiplier  field  of  view,  is  about  twice  as  long  as  indicated  by  the 
photomultiplier.  The  test  time,  as  observed  at  a  given  shock  tube  station 
is  defined  as  the  time  interval  from  the  arrival  of  the  shock  front  to  the 
arrival  of  the  interface.  It  is  also  interesting  to  observe  that  the  ion 
probe  response  reaches  a  steady  level  sooner  than  the  photomultiplier,  and 
then  remains  constant,  despite  the  fact  that  chemical  relaxation  is  still 
proceeding.  The  photomultiplier  field  of  view  is  directly  along  a  shock 
tube  diameter.  Collimation  is  achieved  by  two  0.15  mm  wide  slots  placed 
12  inches  apart.  The  aperture  slots  are  about  3.25  mm  long. 
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In  contrast  to  the  photomultiplier  response  corresponding  to  shock 
conditions  at  1  mmHg  initial  pressures,  the  photomultiplier  trace  at 
0.1  mm  initial  pressure,  shown  in  Fig.  B- 6 ,  does  not  indicate  any  in¬ 
creased  radiation  interface  activity.  Consequently,  at  low  pressure, 
photomultipliers  do  not  give  a  good  indication  of  test  time. 
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FIG.  B-6  TYPICAL  PHOTOMULTIPLIER  RESPONSE 

FOR  0.1  mmHg  INITIAL  SHOCK  TUBE  PRESSURE 
AND  LOW  SHOCK  VELOCITY  AT  THAT  PRESSURE 


Figure  B-7  shows  some  image  converter  photographs.  In  the  upper 
picture,  taken  at  1  mmHg  initial  shock  tube  pressure,  the  shock  interface 
separation  at  the  time  was  always  too  large  to  be  recorded  on  the  same 
shot.  The  picture  indicates  the  flow  field  structure  about  the  flat  end 
of  a  one  inch  diameter  cylinder.  In  the  lower  photograph  at  0.1  mmHg 
initial  pressure,  the  ragged  interface  is  vividly  apparent,  indicating 
almost  no  test  time.  An  ion  probe  sitting  in  the  lower  part  of  the  tube 
would  indicate  test  time  of  the  order  of  10  /xsec. 

Shock  velocity  versus  capacitor  energy  input  into  the  driver  is  shown 
in  Fig.  B-8.  To  date  the  shock  tube  has  been  tested  at  1  mmHg  and 
0.1  mmHg  initial  pressures.  The  capacitor  bank  voltage  variation  ranged 
from  7  kv  to  15  kv.  It  is  perhaps  worth  mentioning  that  our  particular 
driver  triggering  system  has  permitted  us  to  break  down  the  driver  gas  at 
very  low  voltages,  say  6.5  kv ,  and  200  psi  driver  pressures  without 
difficulties. 

Figures  B-9  and  B-10  present  the  observed  shock  tube  test  times  about 
19  ft  from  the  exit  of  the  diverging  section  for  1  mmHg  and  0.1  mmHg 
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FLOW  ABOUT  I  INCH  DIAMETER  CYLINDER, 
FLAT  END  OF  CYLINDER  IS  FACING  FLOW. 
SHOCK  FRONT  NOT  VERY  BRIGHT,  BECAUSE 
U.V.  RADIATION  IS  ABSORBED  IN  PYREX 
GLASS,  BRIGHT  LINE  BEHIND  SHOCK 
FRONT  IS  1/4  -  INCH  BRASS  SUPPORT  ROD 

STATION  719 
P,  =  I  mmHg 
Ut  =  4.5  mm//isec 


NOTE  RUGGED  INTERFACE:  LOCAL  TEST 
TIME  VARIES  FROM  A  ZpfC  MINIMUM  TO 
ABOUT  8/isec  MAXIMUM 

STATION  719 
P,  =  O.l  mmHg 
U#  s  6.5  mm/^sec 

SHOCK  FRONT  HAS  BEEN  ENHANCED 


* 


FIG.  B-7  IMAGE  CONVERTER  PICTURES 
(Courtesy,  Electro  Optical  Corp.) 


SHOCK  VELOCITY,  Us - mm/^i.sec 


FIG.  B-8  SHOCK  VELOCITY  vs.  CAPACITOR  VOLTAGE 

FOR  PRESSURE-DRIVEN,  ARC-HEATED  SHOCK  TUBE 
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TEST  TIME -  (ifc 


SHOCK  VELOCITY,  U,  - nun//iuc 

T»-««»7 -it 


FIG.  B-9  SHOCK  TUBE  TEST  TIMES  vs.  SHOCK  VELOCITY 
AT  1  mmHg  INITIAL  SHOCK  TUBE  PRESSURE 
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FIG.  B-10  SHOCK  TUBE  TEST  TIMES  vs.  SHOCK  VELOCITY 
AT  0.1  mmHg  INITIAL  SHOCK  TUBE  PRESSURE 
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initial  shock  tube  pressures,  respectively.  The  test  times  were  obtained 
(rom  observations  of  photomultiplier  response  characteristics  and  from 
ion  probes.  Photomultiplier  test  times  correspond  to  the  time  between 
the  center  of  the  initial  spiked  deflection  and  the  secondary  increase  in 
radiation  level,  generally  associated  with  the  interface.  Test  times 
obtained  from  ion  probes  were  similarly  counted  from  the  initial  rise  due 
to  the  shock  front  to  the  subsequent  increase  in  activity  when  the  inter¬ 
face  arrives.  At  low  shock  velocities,  however,  no  significant  increase 
in  instrument  deflection  is  observed  by  either  the  photomultiplier  or  the 
ion  probes.  In  those  cases,  an  arbitrary  minimum  fluctuation  level  of 
about  50  percent  is  set;  when  probe  current  fluctuations  exceed  this 
level,  the  interface  is  presumed  to  have  arrived. 

Theoretical  test  time  expectations,  using  Mi  re  1 9 s 23  equations ,  have 
been  superimposed  on  the  experimental  results  in  Figs.  B-9  and  B- 10  for 

an  L/D  ratio  of  19.  The  wide  scatter  of  experimental  data  points  is 

characteristic  of  arc-driven  shock  tubes  (see  for  example,  reports  from 
Avco22  and  JPL24) .  Note  that  a  few  data  points  fall  above  the  theoretical 
curve;  this  is  probably  due  to  the  fact  that  we  based  our  L/D  ratio  on 
the  distance  from  the  diverging  channel  exit  rather  than  on  the  distance 
from  the  diaphragm.  The  diaphragm  is  about  5  ft  upstream  of  the  diverging 
channel  exit,  but  over  this  distance  the  channel  tube  cross  section  varies 
from  a  VA  inch  diameter  to  the  12  inch  test  tube  diameter. 

Since  the  original  shock  tube  length  was  short,  L/D  f  25,  even  the 
calculated  theoretical  test  times  are  short.  To  remedy  the  situation, 
the  L/D  was  recently  increased  to  36.  A  few  test  times  observed  at  the 

new  far  end  position  are  also  plotted  on  the  graph.  Often,  the  photo¬ 

multipliers  indicated  substantially  shorter  test  times  than  did  the 
free-stream  ion  probes  placed  about  1  to  2  inches  into  the  shock  tube. 

This  indicated  a  curved  interface  with  the  flow  in  the  center  of  the  tube 
projecting  forward  by  about  3  to  10  inches  relative  to  the  flow  near  the 
wall.  The  existence  of  a  curved  interface  was  visually  verified  from  the 
image  converter  (see  Fig.  B-7). 

3.  TEST  FOR  EQUILIBRIUM 

Measurements  were  made  at  0.1  an  1  1.0  mmHg  to  check  that  the  slug  of 
ionized  gas  was  indeed  in  ionization  equilibrium.  This  was  done  by  phase 
measurements  made  at  33  Gc  with  a  microwave  interferometer.  The  trans¬ 
mission  path  was  between  a  pair  of  horns  spaced  2  inches  apart  and 
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located  symmetrically  about  the  tube  axis  (see  Fig.  B-ll).  The  horns 
\  were  mounted  so  that  their  apertures  were  flush  with  the  surface  of  the 
plates,  2  inches  from  the  plate  leading  edge.  The  plates  were  2  inches 
wide  and  3  inches  long  with  sharp  tapers  at  their  leading  edges  to  mini¬ 
mize  flow  disturbances.  The  horn  apertures  were  covered  with  Mylar  tape. 

Using  this  system,  electron  densities  were  checked  over  the  range  of 
1010  to  3  x  1012  elec/cc  at  1  mmHg.  Typical  instrument  responses  are 
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FIG.  B-ll  EXPERIMENTAL  CONFIGURATION  OF  MICROWAVE  INTERFEROMETER 
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0.02  v/4m,  1000ft 
v=0.04»ln  ^ 

STATION  755 
R  =  0.1  mmHg 
Ut  s  4.5  mm/u**c 
(n+).  =1.1  a  I01*  iont/cc 
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FIG.  B-12  EQUILIBRIUM  SLUG  ELECTRON  DENSITY  MEASURED  WITH  MICROWAVE 
INTERFEROMETER  AT  0.1  mmHg  INITIAL  SHOCK  TUBE  PRESSURE 
AND  RESPONSE  OF  0.25  mm  dia  P^EE-STREAM  ION  PROBE 
(Upper  Trace:  Interferometer)  * 

(Lower  Trace:  Ion  Probe) 


presented  in  Figs.  B-12  and  B-13.  The  electron  densities  inferred  from 
the  interferometer  measurements  agreed  to  within  30  percent  of  the  values 
of  the  equilibrium  electron  densities  calculated  from  the  measured  data 
on  shock  velocity.  These  results  were  obtained  with  the  tube  fired  over 
thirty  times  without  cleaning. 
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FIG.  B-13  EQUILIBRIUM  SLUG  ELECTRON  DENSITY  MEASURED  WITH  MICROWAVE 
INTERFEROMETER  AT  1  mmHg  INITIAL  SHOCK  TUBE  PRESSURE 
AND  RESPONSE  OF  0.25  mm  dia  FREE-STREAM  ION  PROBE 

(Upper  Trace:  Interferometer) 

(Lower  Trace:  Free-Stream  Ion  Probe) 
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The  tube  was  also  operated  with  dry  air  as  well  as  with  ordinary  room 
air.  No  significant  difference  was  detectable  either  with  probes,  or  with 
the  interferometer.  Neither  was  a  change  in  agreement  with  equilibrium 
calculatio^  detected. 

Measurements  made  at  0.1  mmHg  showed  more  sensitivity  to  tube  clean¬ 
liness  when  fired  at  slower  velocities.  Measurements  made  with  40  shots 
on  the  tube  without  cleaning  indicated  electron  densities  as  much  as  an 
order  of  magnitude  lower  than  equilibrium.  After  cleaning,  the  interfer¬ 
ometer  still  read  lower  electron  densities  than  equilibrium,  but  the  level 
was  raised  to  within  a  factor  of  two  of  equilibrium.  A  0.25  x  6.4  mm 
probe  mounted  between  the  interferometer  plates  indicated  electron  densi¬ 
ties  close  to  those  deduced  from  the  interferometer.  A^similar  probe 
mounted  forward  of  the  interferometer  plates  indicated  electron  densities 
higher  than  the  probe  mounted  between  the  plates.  The  forward  probe 
measured  values  closer  to  equilibrium  than  the  probe  mounted  between  the 
plates. 

These  results  indicate  that  the  plates  were  disturbing  the  flow  at 
0.1  mmHg  initial  pressures  and  causing  a  lower  electron  density  to  appear 
between  the  plates  than  exists  in  free  stream.  The  free-stream  electron 
density  as  measured  by  the  forward  probe  indicated  electron  densities 
within  a  factor  of  two  of  equilibrium.  The  agreement  between  the  inter¬ 
ferometer  and  the  probe  located  between  the  plates  is  taken  as  evidence 
that  the  0.25  mm  free-stream  probe  is  accurately  measuring  the  electron 
density. 

Measurements  made  with  identical  probe  space  120  degrees  apart  around 
the  circumference  of  the  tube  indicate  that  variations  in  electron  density 
around  the  tube  can  be  as  high  as  a  factor  of  two. 
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TEST  FACILITIES  USED  IN  THE  RVIP  EXPERIMENTAL  PROGRAM 
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APPENDIX  C 


TEST  FACILITIES  USED  IN  THE  RVIP  EXPERIMENTAL  PROGRAM 

1.  ELECTROMAGNETIC  SHOCK  TUBE 

The  electromagnetic  shock  tube  used  in  this  investigation  is  shown 
in  Figs.  C- 1  and  C- 2 .  The  discharge  section  consists  of  a  tapered  glass 
tube  with  a  solid  electrode  at  the  small  end  and  a  ring  electrode  at  the 
large  end  where  the  cylindrical  pipe  begins.  The  shock  tube  proper  is  a 
15  cm-ID  Pyrex  pipe  as  shown  in  Fig.  C-2.  The  30  fJ. fd,  20-kv  capacitors, 
the  charging  power  supply,  and  appropriate  switch  gear  are  housed  within 
a  metal  enclosure  (inset).  The  enclosure  also  houses  the  discharge  sec¬ 
tion  as  indicated  in  the  inset  of  Fig.  C-2.  The  heavy  top  plate  to  which 
the  pump  ring  and  cylindrical  test  section  are  attached  also  serves  as 
the  ring  electrode  of  the  discharge  section.  The  switch  between  the 
capacitors  and  the  central  electrode  is  activated  by  reducing  the  pres¬ 
sure  in  the  switch  section  [Fig.  C- l]  until  breakdown  occurs.  The  shock 
tube  is  capable  of  producing  plasma  flows  greater  than  1  cm/usec  and 
electron  densities  greater  than  1014  electrons/ cc .  A  description  of  the 
shock  tube  performance,  including  an  evaluation  of  Langmuir  probe  per¬ 
formance  in  this  environment,  can  be  found  elsewhere.25*26 

For  these  measurements,  the  tube  was  evacuated  to  O.lmmHgand  fired 
at  9  kv.  Under  these  conditions,  the  shock  velocity  is  less  than 
2.0  mm//xsec  and  the  electron  temperature  behind  the  shock  is  estimated 
(from  swept  Langmuir  probe  measurements)  to  be  5000  K.  Under  these  con¬ 
ditions,  the  ion  drift  velocity  (approximately  equal  to  the  shock  ve¬ 
locity)  is  less  than  the  thermal  velocity  of  the  ions  (approximately 
3.5  mm/sec).  Some  of  the  measurements  were  performed  in  the  Region 
milliseconds  after  the  shock  front  passed  so  that  the  ion  drift  velocity 
approached  zero  and  the  electron  temperature  became  300  K  (inferred  from 
swept  Langmuir  probes). 

2.  LOW-PRESSURE  FLAME  FACILITY 

The  low-pressure  chamber  is  a  20-inch  cylindrical  section  of  Pyrex 
glass.  Two  4-inch  diameter  ports  have  been  added  on  opposite  sides  to 
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FIG.  C-2  PHOTOGRAPH  OF  ELECTROMAGNETIC  SHOCK  TUBE 
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provide  flat  glass  viewing  areas.  All  access  to  the  charn 
the  base  plate.  The  steel  plate  forming  the  top  of  the  c 
as  a  short  section  of  the  4-inch  exhaust  line  above  it,  a 
A  water-cooled  baffle  section  to  further  cool  the  gas  was 
exhaust  line  before  the  4-inch  gate  valve  used  for  flow  c 
mixed  e thy lene- oxygen  flame  is  used  in  this  measurement, 
shows  the  test  facility,  including  the  visible  flame.  A 
description  of  the  facility  is  in  Ref.  23. 

For  these  measurements,  the  flame  was  burned  at  7  mnr 
with  potassium  or  sodium  chloride.  The  flame  temperature 
determined  by  a  sodium  line  reversal  measurement.  Microw 
static  probe  measurements  of  the  charged  particle  density 
the  charged  particle  density  is  in  Saha  equilibrium  with 
density  and  temperature.  These  measurements  further  indi 
tence  of  a  negligible  number  of  negative  ions. 
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FIG.  C-3  ETHYLENE-OXYGEN  FLAME  BURNING  AT  8  mmHg  PRESSURE 
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